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Abstract. The present paper commences the study of higher order differ- 
ential equations in composition form. Specifically, we consider the equation 
Lu = div B* V(a div AVtt) = 0, where A and B are elliptic matrices with 
complex-valued bounded measurable coefficients and a is an accretive func- 
tion. Elliptic operators of this type naturally arise, for instance, via a pull- 
back of the bilaplacian A 2 from a Lipschitz domain to the upper half-space. 
More generally, this form is preserved under a Lipschitz change of variables, 
contrary to the case of divergence-form fourth order differential equations. We 
establish well-posedness of the Dirichlet problem for the equation Lu = 0, with 
boundary data in L 2 , and with optimal estimates in terms of nontangential 
maximal functions and square functions. 



1. Introduction 

The last few decades have witnessed a surge of activity on boundary-value 
problems on Lipschitz domains for second-order divergence-form elliptic equations 
— divAVu = 0. Their investigation has, in particular, been guided by two princi- 
ples. 

First, divergence-form equations are naturally associated to a bilinear "energy" 
form, and admit a variational formulation. It turns out that some smoothness of 
the coefficients A in a selected direction is necessary for well-posedness of the un- 
derlying boundary problems in = {(x,t) : x e K", t > 0} (see [14]). This 
observation led to the study of the coefficients constant along a single coordinate 
(the ^-coordinate when n > 2). The well-posedness of the corresponding boundary- 
value problems was established for real symmetric matrices in [55J H2] , and the real 
non-symmetric case was recently treated in [41] [45 , 55 , 36 . In addition, the reso- 
lution of the Kato problem [8] provided well-posedness for complex i-independent 
matrices in a block form; see [6] [47]. Furthermore, a number of perturbation- type 
results have been obtained, pertaining to the coefficient matrices close to "good" 
ones in the sense of the L°° norm [55] [HI El El H2], or in the sense of Carleson 
measures [H [Ml 1301 W\ 113 ESI (Ml E£2 151 El I3H] • 

A seemingly different point of view emerges from the ultimate goal of treating 
boundary-value problems on non-smooth domains rather than just the upper half- 
space. However, it brings to focus equations of the same type as above. Indeed, the 
direct pull-back of the Laplacian A from a Lipschitz domain {(x,t) : t > ip(x)} to 
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M™ yields a boundary- value problem for an operator of the form — div AX7u = 
on K™ +1 , and the corresponding matrix A is, once again, independent of t. More 
generally, if p is a change of variables, then there is a real symmetric matrix A = 
A{x) such that if Am = in f2 and u = u o p, then div AXJu = in 

The model higher order differential operator is the bilaplacian A 2 = — A(— A). 
Investigating the behavior of biharmonic functions under changes of variables, we 
find that there exists a scalar-valued function a and a real elliptic matrix A such 
that, if A 2 u = in SI, then 

(1.1) divAV(a divAVu) = in p _1 (0), 

More generally, such a form is preserved under changes of variables. We em- 
phasize that this is not the case for higher order operators in divergence form 

In addition, the form appearing in (jl.lj) mimics the structure of the bilaplacian 
as a composition of two Laplace operators. As it turns out, this is an important 
feature that underpins several key properties of the solutions to the biharmonic and 
polyharmonic equations (— A) m u = 0, m > 1. 

Motivated by these considerations, the present paper commences the study of 
well-posedness problems for higher order equations in composition form. Specifi- 
cally, consider the equation 

(1.2) div B* V(a div AVu) = 0. 

Here a : R n+1 H > C is a scalar-valued accretive function and A and B are (n + 1) 
x (n + 1) elliptic matrices with complex coefficients. That is, there exist constants 
A > A > such that, if M = A or M = B, then 

(1.3) A < Rea(X) < \a(X)\ < A, A|t?| 2 < Refj ■ M(X)r), |f • M(X)r)\ < A\t]\\£\ 

for all X £ R n+1 and all vectors £, -q £ C™ +1 . The second-order operators div^4V 
and divi3*V are meant in the weak sense; see IDcfinition 2.9l below for a precise 
definition. We assume that the coefficients a, A and B are t- independent; no 
additional regularity assumptions are imposed. For technical reasons we also require 
that A and B satisfy the De Giorgi-Nash-Moser condition, that is, that solutions u 
to divMV-u = are locally Holder continuous for M = A, B, A* and B*. 

The main result of this paper is as follows. We show that whenever the second 
order regularity boundary- value problems for A and B are well-posed, and whenever 
the operator L = divi3*Va divAV is close to being self-adjoint, the L 2 -Dirichlet 
problem 

{divB*V(adivylVu)=0 in 
u = f oncM + +1 , V/£i 2 (M"), 
e n+1 -AX7u = g on 5M" +1 , g £ L 2 (M"), 
has a unique solution u that satisfies the optimal estimates 

N+(X7u) £ L 2 (R n ) and \\\t div AX7u\\\ + < oo, 

where 

1/2 

N+(Vu)(x) =supiY/ |V?f) : y £ R n ,\x - y\ < s) 

l\JB((y,s),s/2) J J 
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and where 

111**1+= / / \F(x,t)\ 2 tdtdx. 

JR™ JO 

Specifically, we will construct solutions whenever ||Ima||Loo( K n) and \\A — -B||i°°(B»i) 
are sufhciently small. It is assumed that the second-order regularity problem 

div AVu = in R£ +1 , u = / on dW^ 1 , N±(Vu) G L 2 (K") 

has a unique solution in both the upper and lower half-spaces whenever V/ € 
L 2 (R"), and that the same is true for A* . By perturbation results in [7J, an analo- 
gous statement is then automatically valid for B and B*. 

We will construct solutions using layer potentials; the De Giorgi-Nash-Moser 
requirement mentioned above is necessary for this approach and at the moment is 
common in the theory of second-order problems. If, for example, A and B are real 
symmetric or complex and constant, then the De Giorgi-Nash-Moser condition is 
valid and regularity problems are well-posed; hence, in these cases there is well- 
posedness of (|1.4|) . We mention that in passing we prove that if the second-order 
regularity problems are well-posed then the solutions necessarily can be written 
as layer potentials^ This fact is new and interesting on its own right. We will 
precisely state the main theorems in ISection 2.21 after the notation of this paper 
has been established. 

Let us point out that to the best of the authors' knowledge, this is the first result 
regarding well-posedness of higher order boundary- value problems with non-smooth 
variable coefficients and with boundary data in LP . For divergence form equations, 
some results for boundary data in Besov and Sobolev spaces LP a , < a < 1, are 
available (see [5J SB]), but they do not reach out to the "end-point" case of LP 
data. Until now, well-posedness in Lipschitz domains with LP data was known 
only for constant coefficient higher order operators (see [30] [53] 54, 64. 58, 591 160] ). 
As explained above, our results extend to Lipschitz domains automatically via a 
change of variables. (See lTheorem 2.341 for the precise statement.) 

Let us discuss the history of the subject and our methods in more detail. We will 
concentrate on higher order operators and only mention the second order results 
that directly affect our methods. The basic boundary- value problem for elliptic 
differential equations of order greater than 2 is the L p -Dirichlet problem for the 
biharmonic operator A 2 . It is said to be well-posed in a domain if, for every 
/ € Wf (dCl) and g £ L p (dfl), there exists a unique function u that satisfies 

(1.5) A 2 w = 0inf7, w = / on 30, v ■ Vu = g on dQ, ]V S! (Vii) e I p (ffl) 
where v is the outward unit normal derivative and 

(1.6) N a (Vu)(X) = sup{|Vtt(F)| : Y e fl, \X - Y\ < (1 + a) dist(Y, dCl)} 

for some constant a > 0. In [56, 16, 62] . well-posedness of the L p -Dirichlet problem 
for A 2 was established in C 1 domains for any 1 < p < oo. (This result is also valid 
in convex domains; see [BHISS].) For general Lipschitz domains, the L 2 -Dirichlet 
problem for A 2 was shown to be well-posed by Dahlberg, Kenig and Verchota in [20] 



This result is tantamount to proving invertibility of the single layer potential. The method of 
establishing injectivity from the regularity problem and from jump relations is known; the authors 
would like to thank Carlos Kenig for bringing this argument to their attention. The method of 
establishing surjcctivity is new and again uses jump relations. 
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(when the domain is bounded; cf. |53[ Theorem 3.7] for domains above Lipschitz 
graphs). 

The sharp range of p for which the L p -Dirichlet problem is well-posed in n- 
dimensional Lipschitz domains is a difficult problem, still open in higher dimensions 
even for the bilaplacian (cf. [TSJ [40] ) . We do not tackle the well-posedness in L p , 
p =/= 2, in the present paper; it is a subject for future investigation. However, let 
us mention in passing that for the bilaplacian the sharp results are only known in 
dimensions less than or equal to 7 [58, 59, 6U] and, in a dramatic contrast with 
the case of the second order boundary- value problems, there is a sharp dimension- 
dependent upper bound on the range of well-posedness. That is, if C R" is a 
Lipschitz domain, then solutions to ()1.5|) are guaranteed to exist only for 2 < p < p n 
for some p n < oo. Related counterexamples can be found in [53[ Theorem 10.7]. 
See also [13] for a review of this and related matters. 

Our methods in the present paper depart from the ideas in [20) and [53j . The 
solution is represented via 

u(X) = -V A f{X) - S A g{X) + £ B ,a,Ah(X). 

Here V A and S A are the classic double and single layer potentials associated to the 
operator div AW , given by the formulas 

V A f(x,t) = - [ e n+1 -A*(y)VTf; t) (y,0)f(y)dy, 
S A g(x,t)= f T£ t) (y,0)g(y)dy 

JR" 

where T^- is the fundamental solution to — divA*V (that is, the solution to 
— divA*Vr^ = Sx)- On the other hand, £b,o,,a is a new layer potential, specifi- 
cally built for the problem at hand, to satisfy 

-a div AW£ B , a ,Ah = l R ^+id 2 l+1 S B 'h. 

See ISection 2.41 This resembles the formula used in 53_ to construct solutions 
to A 2 u = 0. To prove existence of solutions to (|1.4I) or (|1.5I) . in addition to the 
second order results, we require appropriate nontangential maximal function (and 
square-function) estimates for the new potential £s,a,A, as well as the invertibility 
of h i y d n+ i£B } a,Ah\g R rz+i in L 2 . However, beyond the representation formula and 

invertibility argument, our method is necessarily different from |20) and [53j . 

After a certain integration by parts, the bounds on the nontangential maximal 
function of the new potential £h in the case of the bilaplacian become an automatic 
consequence of the Calderon-Zygmund theory and boundedness of the Cauchy in- 
tegral in L 2 . On the other hand, for a general composition operator, the related 
singular integral operators do not fall under the scope of the Calderon-Zygmund 
theory and, because of the presence of non-smooth matrices A and B, are not 
amenable to a similar integration by parts. We develop an alternative argument, 
appealing to some elements of the method in [3], to obtain square function bounds, 
and then employ the jump relations and intricate interplay between solutions in the 
upper and lower half-spaces to obtain the desired nontangential maximal function 
estimates. 

It is interesting to observe that, given an involved composition form of the oper- 
ator, with several "layers" of non-smooth coefficients, the difficulties also manifest 
themselves in the absence of a classical variational formulation. In particular, such 
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standard properties of solutions as the Caccioppoli inequality have to be reproven 
and even the existence of the Green function or fundamental solution in R™ +1 is 
not obvious, in any function space. In the same vein, the existence of the normal 
derivative of a solution cannot be viewed as the result of an integration by parts and 
an approximation scheme. Instead, it once again calls for some special properties 
of the associated higher order potentials. 

Needless to say, our results build extensively on the developments from the theory 
of second-order divergence form operators La — — div AV. We refer the reader to 
[4Ti] for a detailed summary of the theory as it stood in the mid-1990s, and to the 
papers [Ml [44l [24l [45l [5^ E [3 [47j [25l El 13 for more recent developments. 

Finally, let us mention that aside from the the Dirichlet case, it is natural to con- 
sider the Neumann and regularity problems as well as the inhomogeneous equation 
Lu = f, for / in a suitable function space. Recent achievements in this direction 
for higher order equations include [171 [HI] , [S31 [EH EH] , and [TJ [SOI [21 HE] respec- 
tively. Unfortunately, much as in the homogeneous Dirichlet case, they concentrate 
mostly on constant coefficients, with the exception of [2] 148] : these two papers con- 
sider the inhomogeneous problem Lu — f but require that the boundary data have 
extra smoothness in L p . 

The outline of this paper is as follows. In ISection 2l we will define the no- 
tation used throughout this paper and state our main results. In ISection "3l we 
will review known results from the theory of second-order operators of the form 
La = — div AV. In ISection 4[ we will prove fourth-order analogues to some basic 
theorems concerning solutions to second-order equations, such as the Caccioppoli 
inequality. We will construct solutions to (|1.4[) using potential operators and estab- 
lish that these potentials are well-defined and bounded in Sections [5j [HI and [7J The 
invertibility and uniqueness results will be presented in ISection 81 together with the 
end of proof of the main theorems. 



2. Notation and the main theorems 

In this section we define the notation used throughout this paper; in ISection 2.21 
we will state our main theorems. (The proofs will be delayed until ISection 81 ) 

We work in the upper half-space R™ +1 = R n x (0, oo) and the lower half-space 
= R™ x (-oo, 0). We identify <9R± +1 with R". The coordinate vector e = e n+1 
is the inward unit normal to R™ +1 and the outward unit normal to R™ +1 . We will 
reserve the letter t to denote the (n + l)st coordinate in M ,l+1 . 

If O is an open set (contained in R™ or R™ +1 ), then Wf (f2) denotes the Sobolev 
space of functions / e L 2 (VL) whose weak gradient V/ also lies in L 2 (S1), and 
W^^fl) denotes its dual space. The local Sobolev space W^; oc (f2) denotes the 
set of all functions / that lie in Wf (V) for all open sets V compactly contained 
in il. We let W?(R n ) be the completion of Cg°(R n ) under the norm ||/||^2 (En) = 

l|V/|| L 2 (K „ ); equivalently W?(R n ) is the space of functions / G W 1 2 jJoc (R n ) for 
which || V/H/^fRr.) is finite. Observe that functions in W^(M. n ) are only defined up 
to additive constants. 

If u € Wl loc (SY) for some fi C R n+1 , we let Vmu denote the gradient of u in the 
first n variables, that is, V\\u = (dxu, c^u, . . . , d n u). We will occasionally use V|| 
to denote the full gradient of a function defined on R™. 
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As in [3l |6] and other papers, we will let the triple-bar norm denote the L 2 norm 
with respect to the measure (1/|£|) dxdt. That is, we will write 



(2.1) 



|F||| 2 ± = f°° f \F(x, ±t)\ 2 dx \ dt, \\\F\\\ 2 = \\\F\\\l + \\\F\\ 

JO JR™ 1 



with the understanding that a t inside a triple-bar norm denotes the (n + l)st 
coordinate, that is, 



\tF\\\% = [ [ \F(x,±t)\ 2 dxtdt. 

JO JM™ 



We let B(X,r) denote balls in R" +1 and let A(x,r) denote "surface balls" on 
dWp- 1 , that is, balls in R". If Q C E™ or Q C R' i+1 is a cube, we let t{Q) denote its 
side- length, and let rQ denote the concentric cube with side-length r£(Q). If E is a 
set and /i is a measure, we let f denote the average integral j E f dfj, — J E f dfi. 

We will use the standard nontangential maximal function N, as well as the 
modified nontangential maximal function N introduced in |42j . These functions 
are defined as follows. If a > is a constant and x £ M™, then the nontangential 
cone 7±(x) is given by 

(2.2) 1± (x) = {(y,s)eRl +1 :\x-y\<a\s\}. 

The nontangential maximal function and modified nontangential maximal func- 
tion are given by 

(2.3) N ± F(x) = sup{|F(y,s)| : (y,s) G 7±(x)} , 

1 /2 

(2.4) N±F(x) = sup( f / \F\ 2 ) :(y,s)ej±(x) 

I \JB((y,s),\s\/2) J 

We remark that by [27l Section 7, Lemma 1], if we let 

N a F(x) = sup{|^(y,s)| : \x - y\ < as, < s} , 

then for each 1 < p < oo and for each < b < a, there is a constant C depending 
only on p, a and b such that ||iV a F||ip(K«) < C\\NbF\\L P {R n y Thus, for our purposes, 
the exact value of a in (|2.2|) is irrelevant provided a > 0. 

Suppose that A : R™ +1 i— > c("+ 1 ) x ("+ 1 ) is a bounded measurable matrix- valued 
function defined on R™ +1 . We let A T denote the transpose matrix and let A* denote 
the adjoint matrix A T . Recall from the introduction that A is elliptic if there exist 
constants A > A > such that 

(2.5) X\n\ 2 ^Refj-AiXH |£ • A(X)rj\ < A\ v \\t\ 

for all X € R" +1 and all vectors r/, £ £ C" +1 . We refer to A and A as the ellipticity 
constants of A. Recall also that a scalar function a is accretive if 

(2.6) A < Rea(A) < \a(X)\ < A for all X £ R n+1 . 
We say that a function or coefficient matrix A is t-independent if 



(2.7) 



A(x, t) = A(x, s) for all x £R n and all s, t £ R. 
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2.1. Elliptic equations and boundary-value problems. If A is an elliptic ma- 
trix, then for any u £ / oc (fi), the expression Lau = — divAVu £ W 2 1 loc (£l) is 
defined by 

(2.8) (ip,L A u) = I 7pL A u= [\7^ ■ A\7u for all <p £ C™(fl). 

Jn J 

If A and B are elliptic matrices and a is an accretive function, we may define 
L* B (a Lau) in the weak sense as follows. 

Definition 2.9. Suppose u £ W 2 loc (£l). Then Lau = — div AVu is a well-defined 
element of ; oc (^)- Suppose that aL^u = v, for some v £ Wf ; oc (fi), in the 
sense that 



y Vtp • AVu = J Tp^-v for all 93 £ C%°(Q). 



If u satisfies 



J Vv ■ BVn = Jvf for all V € C °°(f7), 

that is, if — div B*Vv = / in the weak sense, then we say that L* B (a Lau) = f. 

Suppose that U is defined in MJ +1 . We define the boundary values of U as as 
the L 2 limit of U up to the boundary, that is, 

(2.10) U\ ml+1 =F provided \jm ± \\U( ■ ,t) - F\\ L 2 {Rn} = 0. 

Given these definitions, we may define the Dirichlet problem for the fourth-order 
operator L* b (ciLa) as follows. 



Definition 2.11. Suppose that there is a constant Co such that, for any / £ 



Wf(R n ) and any g £ L 2 (R n ), there exists a unique function u £ W^^W^ 1 ) that 



satisfies 

r T*L( n T..„.\-n i n R + 



(2.12) 



L* B (aL A u) = hiR™ +1 

u = f, V||U = V/ on9IR r ; +1 , 

e-AVti = g ondlR™ +1 , 

I ||iV + (V u )|| L 2 (R „) + |||tL A «|||+ < C \\Vf\\ L 2 {Rn) + C \\g\\ L ^) 

where L* B (a Lau) — in the sense of Definition 2.91 and where the boundary values 
are in the sense of (|2.10p . 

Then we say that the L 2 -Dirichlet problem for L* b (o,La) is well-posed in R^ +1 . 

We specify V||U = V/ as well as u = / in order to emphasize that u(- ,t) —> f 
in Wl (R n ) and not merely in L 2 (W l ). 

Remark 2.13. The solutions u to (|2.12p constructed in the present paper will also 
satisfy the square-function estimate 

(2.14) |||tV<9HII+ < Cx\\Vf\\l^ n) + CxWgWl^ 

for some constant C\. If u is a solution to a fourth-order elliptic equation with 
constant coefficients, then by [19] we have a square function bound on the com- 
plete Hessian matrix V 2 u. However, in the case of solutions to variable-coefficient 
operators in the composition form of IDcfinition 2.9"! we do not expect all second 
derivatives to be well-behaved, and so (|2.14|) cannot be strengthened. 
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Our main theorem is that, if a, A and B satisfy certain requirements, then the 
fourth-order Dirichlet problem (|2.12j) is well-posed. We now define these require- 
ments. We begin with the De Giorgi-Nash-Moser condition. 

Definition 2.15. We say that a function u is locally Holder continuous in the 
domain ft if there exist constants H and a > such that, whenever B(Xo, 2r) C O, 
we have that 

(2.16) \u(X)-u(X')\ <h' 




_> 



for all X, X' G B(Xo,r). If u is locally Holder continuous in B(X,r) for some 
r > 0, then u also satisfies Moser's "local boundedness" estimate 

(2.17) \u(X)\<c(f \u\ 2 ) 

\JB(X,r) J 

for some constant C depending only on H and the dimension n + 1. 

If A is a matrix, we say that A satisfies the De Giorgi-Nash-Moser condition if 
A is elliptic and, for every open set and every function u such that div AS7u = 
in f2, we have that u is locally Holder continuous in f2, with constants H and a 
depending only on A (not on u or f2). 

Throughout we reserve the letter a for the exponent in the estimate (|2.16p . We 
will show (see Corollary 4.5| below) that if A, A* and B* satisfy the De Giorgi- 



Nash-Moser condition then solutions to L* B (a Lau) = are also locally Holder 
continuous. 

We say that the I/ 2 -regularity problem (R)^ is well-posed in Rl^ 1 if, for each 
/ G Wi(R n ), there is a function u, unique up to additive constants, that satisfies 

divAVu = inR3+ 1 , 

{R)t{ u = f ondR n ± +1 , 

|JV±(Vti)||L*(R»} < C\\Vf\\ L 2 (Rn) . 



Remark 2.18. If N±(Wu) G L 2 (W l ), then averages of u have a weak nontangential 
limit at the boundary, in the sense that there is some function / such that 

\ V 2 

\u(y, s) - f(x*)\ 2 dyds < C\t\N±{Vu)(x*) 

I B((x,t),\t\/2) J 

for all (x,t) G 7±(x*). See the proof of [42, Theorem 3.1a]; here C is a constant 
depending only on the constant a in the definition (|2.2| of j± . If u is locally Holder 
continuous then this implies that u itself has a nontangential limit, which by the 
dominated convergence theorem must equal its limit in the sense of (|2 . 1 0[) . Thus, 
the requirement in (|2. 12[) or (i?)^ 1 that u = f on cM!^ 1 in the sense of vertical 
L 2 limits is equivalent to the requirement that u = f on 9IRJ +1 in the sense of 
pointwise nontangential limits almost everywhere. 

2.2. The main theorems. The main theorems of this monograph are as follows. 

Theorem 2.19. Let a : W l+1 i-> E and A : R n+1 i-> c(" +1 ) x ( n+1 ), where n+1 > 3. 
Assume that 

• a is accretive and t -independent. 
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• a is real-valued. 

• A is elliptic and t -independent. 

• A and A* satisfy the De Giorgi-Nash-Moser condition. 

• The regularity problems (R)^ and (i?)^ 1 are well-posed in RJ +1 . 

Then the L 2 -Dirichlet problem for L* A (aLA) is well-posed inWl +1 , and the con- 
stants Co and C\ in (|2.12[) and (|2.14l) depend only on the dimension n + 1, the 
ellipticity and accretivity constants X, A in (|2.5I) and (|2.6[) . the De Giorgi-Nash- 
Moser constants H and a, and the constants C in the definition of (R)f and (R)£ ■ 

We can generalize IThcorcm 2.191 to the following perturbative version. 

Theorem 2.20. Let A be as in \Theorem 2.1 (A and let a : W L+1 i— > C be accretive 

and t -independent. Let B : W n+1 i-> c(™ +1 ) x (™ +1 ) be t -independent. 

There is some e > 0, depending only on the quantities listed in I Theorem 2.19\ 
such that if 

||Ima||ioc( H n) < £ and \\A - B||i«,( E n) < e, 

then the L 2 -Dirichlet problem for L* B (a La) is well-posed mR+ +1 , and the constants 
Co and C\ in (|2.12p and (|2.14j) depend only on the quantities listed in \ Theorem 2. 1 ffl 

We will see that if e is small enough then B also satisfies the conditions of 
ITheorem 2.19( see IThcorcm 2.2l"1 It is possible to generalize from R™ +1 to domains 
above Lipschitz graphs; see IThcorcm 2.341 below. 

We remark that throughout this paper, we will let C denote a positive constant 
whose value may change from line to line, but which in general depends only on 
the quantities listed in IThcorcm 2.191 any other dependencies will be indicated 
explicitly. 

In the remainder of this section we will remind the reader of some known sufficient 
conditions for a matrix A to satisfy the De Giorgi-Nash-Moser condition or for the 
regularity problem (R)f to be well-posed. To prove Theorems 12.191 and 12.201 we 
will need some consequences of these conditions; we will establish notation for these 
consequences in ISection 2.31 

We begin with the De Giorgi-Nash-Moser condition. Suppose that A is elliptic. 
It is well known that if A is constant then solutions to divAVw = are smooth 
(and in particular are Holder continuous). More generally, the De Giorgi-Nash- 
Moser condition was proven to hold for real symmetric coefficients A by De Giorgi 
and Nash in [23( [52] and extended to real nonsymmetric coefficients by Morrey in 
[51] , The De Giorgi-Nash-Moser condition is also valid if A is t-independent and 
the ambient dimension n + 1 = 3; this was proven in [3l Appendix B]. 

Furthermore, this condition is stable under perturbation. That is, let Aq be ellip- 
tic, and suppose that Aq and Aq both satisfy the De Giorgi-Nash-Moser condition. 
Then there is some constant e > 0, depending only on the dimension n + 1 and 
the constants A, A in (|2.5|) and H, a in ()2 . 16|) . such that if ||j4 — j 4 ||l°°(r™+ 1 ) < £ > 
then A satisfies the De Giorgi-Nash-Moser condition. This result is from [3]; see 
also [TQl Chapter 1, Theorems 6 and 10]. 

We observe that in dimension n+1 > 4, or in dimension n+1 = 3 for i-dependent 
coefficients, the De Giorgi-Nash-Moser condition may fail; see [3TJ for an example. 

The regularity problem (i?)^ 1 has been studied extensively. In particular, if A is 
t-independent, then {R)^ is known to be well-posed in K!^ +1 provided A is constant, 
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real symmetric ([42]), self- adjoint ([7]), or of "block" form 

A(x) - ( A & ° 

for some n x n matrix Ay and some complex- valued function a±. The block case 
follows from validity of the Kato conjecture, as explained in 10. Remark 2.5.6]; 
see [8] for the proof of the Kato conjecture and [11] Consequence 3.8] for the case 
o±#l. 

Furthermore, well-posedness of (i?)^ 1 is stable under perturbation by [Jj; that is, 
if {R)f° is well-posed in RJ +1 for some elliptic ^independent matrix A , then so is 
{R)2 for every elliptic i-independent matrix A with \\A — Aq\\l^ small enough. 

We mention that if A is a nonsymmetric matrix, the L 2 -regularity problem need 
not be well-posed, even in the case where A is real. See the appendix to [45] for a 
counterexample. 

We may summarize the results listed above as follows. 

Theorem 2.21. Let A be elliptic and t -independent, and suppose that the dimen- 
sion n+1 is at least 3. If A satisfies any of the following conditions, then A satisfies 
the conditions o ftTheorem 2.1$ 

• A is constant. 

• A is real symmetric. 

• A is a 3 x 3 self-adjoint matrix. 

• A is a real or 3 x 3 block matrix. 

If Ao satisfies the conditions of I Theorem TBj and if || A- A || i» (r) < e for 
some e depending only on the quantities enumerated in \ Theorem 2.19i then A also 
satisfies the conditions of \Theorem 2.1(K 

2.3. Second-order boundary-value problems and layer potentials. In or- 
der to prove ITheorem 2.19| we will need well-posedness of several second-order 
boundary-value problems and good behavior of layer potentials; these conditions 
follow from well-posedness of (i?)^ 1 and (i?)^ . 

We say that the oblique L 2 -Neumann problem {N^)^ is well-posed in R± +1 if, 
for each g e L 2 (R"), there exists a unique (modulo constants) function u that 
satisfies 



divAVu = inR^ +1 , 



d n+1 u = g on<9RJ +1 , 

\N±{Wu)\\ L 2 (Rn} < C\\g\\ L 2 {R ny 

By [51 Proposition 2.52] and [71 Corollary 3.6] (see also [7J Proposition 4.4]), if A is 
t- independent, then (i?)^ 1 is well-posed in RJ +1 if and only if (N 1 -)^ is well-posed 



We say that the L 2 -Dirichlet problem (D)^ is well-posed in R± +1 if there is 
some constant C such that, for each / G L 2 (R"), there is a unique function u that 
satisfies 

divAVw^O inRJ +1 , 
(D)U u = f oncM'i +1 , 

k \\N ± u\\ LW <C\\f\\ L2(Kn) . 
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Observe that if u is a solution to (N ± )f with boundary data /, then v = d n +iu is 
a solution to (Z?)^ 1 with boundary data /; thus, well-posedness of (N- 1 )^ implies 
existence of solutions to (D)^- Recall that well-posedness of (N^)-^ also implies 
existence of solutions to (i?)^ 1 . It is possible to show that if (i?)^ 1 solutions exist 
then solutions to (D)% are unique; see the proof of [3l Lemma 4.31]. 

Thus, if (R)f is well-posed in K!^ +1 then so is (D)f. We observe that this result 
was proven for real symmetric A by Kenig and Pipher in |42j . 

The L 2 -Neumann problem (N)^ more usually considered differs from {N- 1 )^ in 
that the boundary condition is e • AS/u = g rather than d n +\u = g on cM™ +1 . We 
remark that if A is constant, self-adjoint or of block form then (N)^ is known to be 
well-posed (again see [Jj or [40l Remark 2.5.6] and [8j[Tl]). (N)^ is in many ways 
more natural than the oblique Neumann problem; however, we will not use well- 
posedness of the traditional Neumann problem and so do not provide a definition 
here. 

A classic method for constructing solutions to second-order boundary- value prob- 
lems is the method of layer potentials. We will use the double and single layer 
potentials of the second-order theory, as well as a new fourth-order potential (see 
IScction 2.4[) to construct solutions to (|2.12l) ; thus, we will need some properties of 
these potentials. 

These potentials for second-order operators are defined as follows. The funda- 
mental solution to La = — div A\7 with pole at X is a function such that 
(formally) — div A\7T X — Sx ■ For general complex coefficients A such that A and 
A* satisfy the De Giorgi-Nash-Moser condition, the fundamental solution was con- 
structed by Hofmann and Kim. See [37J Theorem 3.1] (reproduced as lThcorcm 3.121 
below) for a precise definition of the fundamental solution. 

If / and g are functions defined on R™, the classical double and single layer 
potentials T>aI and SAg are defined by the formulas 

(2.22) V A f(x,t) = -f e-A T (y)\7T^ t) (y,0)f(y)dy, 

(2.23) S A g(x,t)= f Tf"(y,0)g(y)dy. 

For well-behaved functions / and g, these integrals converge absolutely for x € K™ 
and for t ^ 0, and satisfy div AW A f = and div^VS^S = in R n+1 \ R n ; see 
IScction 3.21 

We define the boundary layer potentials 

T>aS = ZWlggj+i, (V5 j4 ) ± 5 = VS A g\ m n+i, 

$a9 = ^Ag\ m n+i, Sa g = d n+1 SAg\ m n+i 

where the boundary values are in the sense of (|2.10p . 

We remind the reader of the classic method of layer potentials for construct- 
ing solutions to boundary-value problems. Suppose the nontangential estimate 
\\N±(\7S A g)\\ L 2 {Rn) < C\\g\\ L 2 {Rn) is valid. Then if is invertible L 2 (R n ) h> 
W?(R n ), then u = SA^S^^g) is a solution to (R)^ with boundary data g. 
Similarly, if S^ :± or the operator g i— > e ■ AftSA^g is invertible on L 2 (M n ), 
then we may construct solutions to (N ± )f or (N)%, respectively. The adjoint to 
g i-> e- AiySA^g is the operator , and so if we can construct solutions to (N)f 
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using the single layer potential then we can construct solutions to (D)^ using the 
double layer potential. In the case of i-independent matrices, we may also construct 
solutions to (D)2 by using invertibility of S^' . 

Thus, if layer potentials are bounded and invertible, we have well-posedness of 
boundary-value problems. The formulas ()2.22|) and (|2.23|) for solutions are often 
useful; thus, the layer potential results above are of interest even if well-posedness 
is known by other methods. It turns out that we can derive the layer potential 
results from well-posedness. 

Lemma 2.24. Suppose (R)£ and (R)f are well-posed in RJ +1 . Then there exists 
a constant C such that 

(2.25) \\\td 2 t S A g\f= f \d 2 t S A g{x,t)\ 2 \t\dxdt<C\\g\\ 2 mm . 

Proof. Recall that if (i?)^ 1 is well-posed then so is (D)^ ■ The square-function 
estimate (|2.25j) follows from well-posedness of (D)^ and (i?)^ via a local T(b) 
theorem for square functions. This argument was carried out in [3j Section 8] in 
the case where A is real and symmetric; we refer the reader to [34, Section 5.3] for 
appropriate functions b to use in the general case. (The interested reader should 
note that [Ml Chapter 5] is devoted to a proof of (|2.25[) in the case of elliptic 
systems.) □ 

It was observed in [38l Proposition 1.19] that by results from [3] and [5], if (|2.25[) 
is valid then 

(2.26) \\N±(VS A g)\\ LHMn) < C\\g\\ L 2 {Rn) . 

Using classic techniques involving jump relations, we will show fsee lTheorem 3.271 
below) that if (R)£ is well-posed in E^ +1 then is invertible L 2 (R n ) i-> W?(R n ), 
and if in addition (N ± )2 is well-posed then S^'^ is invertible. Thus, if (i?)^ 
and (R)2 are well-posed, then not only do solutions to (i?)^, (N- 1 )^ and (D)^ 
exist, they are given by the formulas u = SA{{S A z )~ 1 g), u — SA((S A ~' ± ) 1 g) and 
u = 9 n+ i5 > i((5^' ± ) _1 fli), respectively. 

In many of the results in the later parts of this paper, we will only need a few 
specific consequences of well-posedness of (i?)^ 1 and (i?)^ 1 . These consequences are 
as follows. 

Definition 2.27. Suppose that A and A* are elliptic, i-independent and satisfy 
the De Giorgi-Nash-Moser condition. 
Suppose that 

• the square- function estimate (|2.25[) is valid, 

• the operators Sj; ,=t are invertible L 2 (R n ) \-> L 2 (W l ), and 

• solutions to the Dirichlct problem (D)^ are unique in and in R" +1 . 
Then we say that A satisfies the single layer potential requirements. 

We remark that if A satisfies the single layer potential requirements then (N- 1 )^ 
(and hence (R)f and (D)^) are well-posed, and so both A and A* satisfy the single 
layer potential requirements if and only if A satisfies the conditions of lThcorcm 2.191 
We will prove a few bounds under the assumption that A (not necessarily A*) 
satisfies the single layer potential requirements. We also remark that because we 
have no need of well-posedness of the Neumann problem (N)^, we have not required 
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invertibility of T>^ or its adjoint. Consequently, our layer potential requirements 
are weaker than those considered elsewhere in the literature. 

2.4. Layer potentials for fourth-order differential equations. We will con- 
struct solutions to the fourth-order Dirichlet problem (|2.12l) using potential oper- 
ators. For u = £B.a,Ah to be a solution to L* B (a LA{£h)) — 0, we must have that 
v = a LA£s,a,Ah is a solution to Lb*v = in R™ +1 . We choose v = c^ +1 <Sb»/i; 
if B* is i-independent then L* B (d^ +1 SB'h) = d% +1 L B (SB*ti) — 0. It will be seen 
that with this choice of v, the operator Sb.u.a is bounded and invertible in some 
sense. 

Formally, the solution to the equation 

— a div A\7£ B , a ,Ah = l R ^+id^ +1 SB'h 

is given by 

£ B ,a,Ah{x,t) = I I T( ) {x 1 t)-^—d 2 s S B -h{y 1 s)dsdy. 
At- Jo 1 ; a (y) 

However, to avoid certain convergence issues, we will instead dchnc 

(2.28) £ B , a ,Ah(x, t) = F B>a ,Ah(x, t) - S A ^B* + hj (x, t) 
where the auxiliary potential J~s,a,A is given by 

(2.29) F B ,a, A h{x, t) = ~f J d s T (y,s) ( x > *) i d « S B* Kv, s ) ds d V- 

That these two definitions are formally equivalent may be seen by integrating by 
parts in s. 

In ISection 5l we will show that the integral in (|2.29[) converges absolutely for 
sufficiently well-behaved functions h; we will see that if A and B* satisfy the sin- 
gle layer potential requirements, then by the second-order theory the difference 
^B.a.Ah — J-B.a.Ah = — 5^((1 / 'a)S B ^ + ti) is also well-defined and satisfies square- 
function and nontangential bounds. 

2.5. Lipschitz domains. As we discussed in the introduction, by applying the 
change of variables (x, t) i— > (x, t — ip(x)), boundary-value problems for the second- 
order operator div A\7 in the domain 57 given by 

(2.30) Q = {(x, t) : x G R n , t > (f(x)} 

may be transformed to boundary-value problems in the upper half-space for an 
appropriate second-order operator div AV. In particular, the theory of harmonic 
functions in domains of the form (j2.30[) is encompassed by the theory of solutions 
to div AVu = 0, for elliptic t-independent matrices A, in the upper half-space. 

We now investigate the behavior of fourth-order operators under this (or another) 
change of variables. Let p : Q M- be any bilipschitz change of variables and let 

J p be the Jacobean matrix, so V(wo p) = jJ(Vu) o p. For any accretive function a 

and elliptic matrix M, we let a and M be such that 

(2.31) a = 7^f' J p MJj = \J p \(Mop). 
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By the weak definition ()2.8[) of L^j = — div MV and elementary multivariable 
calculus, we have that if is a domain, and if u G W 2 / oc (R+ +1 ) and Lj^u G 

Ll c (R n + +1 ), then 

(2.32) L M (uap) = \J p \(L s u)op inn = p- 1 (Rl +1 ). 

Observe that u is Holder continuous if and only if u is; thus M satisfies the De 
Giorgi-Nash-Moser condition if and only if M does. 

Now, suppose that L*~(aL^u) — in K" +1 , where a, A, B are given by (|2.3ip . 
Let u = u o p. By Definition 2Jl v = aL^u lies in Wy oc (K+ +1 ). By (|2"32l> , 
vo p = a {Lav) G W / i 2 ; oc (^) and L B *{v o p) = 0. Thus a (I^u) € Wij oc (£2) and so 
L* B (a Lau) is well-defined in O, and furthermore L* B (a Lau) — 0. Thus, existence 
of solutions in O follows from existence of solutions in 

Similarly, if L* b (o,Lau) = in and u = u o then L~(aL^{t) = in 
Thus, uniqueness of solutions in O follows from uniqueness of solutions in K/T 1-1 . 

We remark that the preceding argument is valid with O and K™ +1 replaced by 
V and p(V) for any domain V. 

Because we wish to preserve i-independence, we consider only the change of 
variables p(x,t) — [x,t — <p(x)). This change of variables allows us to generalize 
IThcorcm 2.201 to domains O of the form (|2.30|) . The argument is straightforward; 
however, to state the result we must first define the fourth-order Dirichlet and 
second-order regularity problems in such domains. 

Let O be a Lipschitz domain of the form (|2.30[) . Let v denote the unit outward 
normal to O and a denote surface measure on 90. Let W?(dQ.) denote the Sobolev 
space of functions in L 2 (90) whose weak tangential derivative also lies in L 2 (90); 
in both cases we take the norm with respect to surface measure. We say that U = F 
on 90 if F is the vertical limit of U in L 2 (90), that is, if 

(2.33) lim / \U(X + te n+1 ) - F(X)\ 2 da(X) = 0. 

If / is defined on 90, we let V T / be the tangential gradient of / along 90. If u is 
defined in O, then Vnu is the gradient of u parallel to 90; that is, if X = Xq +te n +\ 
for somclo € 90 and some t > 0, then V\\u(X) = Vu(X) - (u(X ) ■ Vu(X))u(X ). 

If O is the domain above a Lipschitz graph, we say that the L 2 -Dirichlet problem 
for L b (clLa) is well-posed in O if there is a constant Cq such that, for every 
/ G W 2 (dn) and every g G L 2 (90), there exists a unique function u that satisfies 

L b (ciLau) = in O, 

u = f, V || it = V T / on 90, 
v ■ AVu = g on 90, 

< 

\\ N n(Vu)\\ L 2 {dn) < Cb||V T /||x,2(R») + C \\g\\ L 2 {Rn) , 

[\L A u(X)\ 2 dist(X,dn)dX <C \\V T f\\ L 2 (Rn) + C \\g\\ L 2 (Rn) . 
v Jn 
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The modified nontangential maximal function Nfi is given by 

1/2 

N n (Vu)(X) = sup} ( f \Vu\ 2 



I \J B(Y,dist(Y,dQ)/2) 

:Yefl, \X - Y I < (1 + a) dist(T, dfl) 

As in lRcmark 2.18| we also have that u — > f nontangcntially in the sense that 
lim u(Y) = f(X), j(X) = {Y G fi, \X-Y\ < (1 + a) dist(Y, dSl)} . 

We say that the L 2 -regularity problem (i?)^ is well-posed in f2 if for every / G 
Wf (dCt) there is a unique function it that satisfies 

!divAVu = in O, 
u = / on 90, 
||^n(V M )|| L 2 (R „ ) < C||V T /|| L 2 (R „) 

where u = f in the sense of either (|2.33|) or in the sense of nontangential limits. 

^ 4 



Clearly, (R)f is well-posed in fi if and only if (R)£ is well-posed in ™> n+1 



We may now generalize IThcorcm 2.20l to Lipschitz domains. 

Theorem 2.34. Let ft = {(x,t) : t > <p(x)} for some Lipschitz function (p. Let 
a : M" +1 h> C, and let A, B : R" +1 h> c( n+1 ) x ( n+1 ), where n + 1 > 3. 

Suppose that a, A and B are t-independent, that a is accretive, that A and A* 
satisfy the De Giorgi-Nash-Moser condition, and that (i?)^ and (i?)^ 1 are well- 
posed in f2 and in f2 . 

Then there is some e > 0, depending only on ||V<^||x,oo™n) and the quantities 
listed in \Theorem 2.1§\ such that if 

Hlmallioc^n) < e and \\A - B\\ L ^ {Mn) < e, 

then the L -Dirichlet problem for L* b {clLa) is well-posed in fi. 

3. Preliminaries: the second-order theory 

In this section, we will review some known results concerning solutions to second- 
order elliptic equations of the form div AVu = 0, and more specifically, concerning 
solutions to second-order boundary- value problems. In Sections 13.11 and 13.21 we 
will discuss the second-order fundamental solution and some properties of layer 
potentials. 

In this section we will state several results valid in M™ +1 ; the obvious analogues 
are also valid in R™ +1 . 

The following two lemmas are well known. 

Lemma 3.1 (The Caccioppoli inequality). Suppose that A is elliptic. Let X G 
W l+1 and let r > 0. Suppose that div AVu = in B(X,2r), for some u G 
W?(B(X,2r)). Then 

ivd 2 < S ! \u? 



I B(X,r) 1 J B(X,2r)\B(X,r) 

for some constant C depending only on the ellipticity constants of A and the di- 
mension n + 1. 
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Lemma 3.2. [491 Theorem 2]. Let A be elliptic, and suppose that divAVu = 
in B(X,2r). Then there exists a p > 2, depending only on the constants X, A in 
(|2.5p . such that 

\Vu\A <c(i \Vu\ 2 \ . 

>B(X,r) J \JB(X.2r) J 

These conditions may be strengthened in the case of t-independent coefficients. 
Suppose that Q C K" isacube. If u satisfies div A\7u = in 2Qx (t-£(Q), t+i{Q)), 
and if A is t- independent, then by O Proposition 2.1], there is some po > 2 such 
that if 1 < p < po, then 

/ r \Vp / r rt+e(Q)/i x 1/2 

(3.3) If \Vu{x,t)\ p dx) <C(f f \Wu(x,s)\ 2 dsdx) . 

\Jq J \J2Q Jt-t(Q)/i J 

We will show that a similar formula holds for solutions to fourth-order equations 
in ILcmma 4.81 below. 

More generally, we have the following lemma. 

Lemma 3.4. Suppose f S L 2 (M. n *-> C" +1 ) is a vector-valued function, and that 
div AS7u — in 2Q x (t — £(Q), t + £{Q)) for some t-independent elliptic matrix A. 
Then 

r _ r rt+t(Q)/2 

(3.5) f \Vu(x,t) - f(x)\ 2 dx <C f f \Vu(x,s)- f(x)\ 2 dsdx. 

JQ J2Q Jt-i(Q)/2 

Proof. Let v(x,t) = f^^qy^ u{x, s) ds. Then 

\Vu(x,t) - f(x)\ 2 dx <2f \Vu(x,t) - Vv{x,t)\ 2 dx+2 f \Vv(x,t) ~ f{x)\ 2 dx. 

Define D(r) = 2Qx(t-r£{Q),t+r£{Q)). If A is t-independent then div AV(u-v) = 
in £)(3/4). Applying f|3 . 3[) to u — v, we have that 



/ |V«(- ,t)-Vv(- ,t)\ 2 <C -f |Vu-Vw| 



v u — v v 1 2 

By definition of v, 

rt+e(Q)/4 



\Vv(x,t) - f(x)\ 2 dx 



Vu(x, s) — f(x) ds 
t-e(Q)/4 



dx. 



Writing Vu — Vi> = (Vw — /) + (Vi> — /), and again applying the definition of v 
to bound Vv(x, s) — f completes the proof. □ 

We also have the following theorems from [5] . Although we quote these theorems 
for t-independent coefficients only, in fact they are valid for t-dependent coefficients 
that satisfy a Carleson-measure condition. 

Theorem 3.6. Theorem 2.4(i)]. Suppose that A is t-independent and elliptic, 
that div,4Vu = in E+ +1 , and that u satisfies the square-function estimate 



[ [ \Vu(x,t)\ 2 tdxdt < oo. 

JQ JR n 
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Then there is a constant c and a function f G L (M n ) such that 

lim||u( •,£)-/- c|| L 2 (R „) =0 
t->o+ 

and such that 

\\N+(u- c)\\ 2 L 2 (Rn) <C f [ \Wu(x,t)\ 2 tdxdt 



where C depends only on the dimension n + 1 and the ellipticity constants X, A 
of A. If A satisfies the De Giorgi-Nash-Moser condition then we may replace N+u 
by N+u. 

In the t-independent setting, (|3.5[) lets us state [3 Theorem 2.3(i)] as follows. 

Theorem 3.7. Suppose that AivAVu = in M" +1 and that N+{Vu) G L 2 (W l ), 
where A is elliptic and t-independent. Then there exists a function G : R™ i— > C" +1 
with ||G||l 2 (K") < C||-^+(Vu)||i2( R ^) such that 

Km\\Vu( ■ ,t) - (5||| 2(Rn) = = Urn ||V«( • , t)f L2{Mn) . 

By the divergence theorem, there is a standard weak formulation of the boundary 
value e • AVu for any solution u to div AVu = with Vw G L 2 (R" +1 ). IThcorem 3~7l 
implies that if u is an appropriately bounded solution, then the boundary value 
e • A\7u\ m ™+i exists in the sense of L 2 limits. By the following result, these two 

formulations yield the same value. 

Theorem 3.8 ([3J Lemma 4.3]). Suppose that A is elliptic, t-independent and 
satisfies the De Giorgi-Nash-Moser condition. Suppose that div^lVw — in 
andN+(Vu) G L 2 (R n ). 

Then there is some function g € L 2 (M. n ) such that, if tp G Cq°(W"" +1 ) ! then 

S7(p(x,t) ■ A(x)\7u(x,t) dx dt = / (p(x,0) g(x) dx. 

n+i y R „ 

Furthermore, — e ■ AVu( ■ ,t) — > g as t -> 0+ in L 2 (K ra ). 

Given a function w that satisfies div AVu — in R!£ +1 , we will frequently analyze 
the vertical derivative d n +\u. In order to make statements about u, given results 
concerning d n+ \u, we will need a uniqueness result. 

Lemma 3.9. Let A be elliptic and t-independent. Suppose that div AVu = and 
d n+ iu = in Rl +1 . 

If there is some constant c > and some to > such that 



(3.10) / 



IVuI 2 < cT 



/or a/Z t > to and all x G W 1 , then u is constant in 10 " 



B((x,t),t/2) 

If N + (Wu) G L 2 (W l ) then (|3~TU1) is valid for all t > 0. Thus in particular, 
uniqueness of solutions to (D)^ implies uniqueness of solutions to (N^)^- 

Proof. Since d n +iu(x,t) = 0, we have that u(x,t) = v(x) for some function v : 
R n i-> C. By letting t -> oo in (pTOjl we have that Vv G L 2 (R"). 
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Choose some x G R™ and some t > 0. If t > Ct, then let Y(t) be the cylinder 
A(x,t) x (r — t/C,t + t/C). By Holder's inequality and lLemma 3.2| there is some 
p > 2 such that 

|Vu| 2 = -/ |Vu| 2 < ^(fr) 1 - 2 ^ / |Vu|P] 



A(x,i) r JY(t) t \J B((x,t),t/C) j 

J B((x,t),2t/C) 

But if t is large enough, then r" r ) 2 r/c) l^ M l 2 — c - By letting r — » oo, we see 
that Vw = and so u{x,t) — v{x) is a constant. □ 

We will also need the following uniqueness result. 

Lemma 3.11. Let A be elliptic, t-independent and satisfy the De Giorgi-Nash- 
Moser condition. Suppose that u + and U- are two functions that satisfy 



div AVu± =0 in RJ +i , N±(Vu±) G L (R™), Vu+\ m n+i = Vu_| aE „+i. 
Then u + and U- are constant in RJ +1 . 

Proof. Let u — u + in R" +1 , w = u_ + c in R™ , where c is such that u + = li- 
on 1™. We claim that div AVu = in the whole space R ,l+1 in the weak sense 
of (HH). This is clearly true in each of the half-spaces R± +1 . Let <p G C,f (R n+1 ). 
Then 



r 



V<^ • AVu = I V<y3 • AVu + I Vip- AVu. 



Let G — Vu\ dRn +i = Vu| aR n+i. Bv lTheorem 3.81 we have that 

W(p(x, t) ■ A(x) Vm(i, t) dx dt = =F / tp{x, 0) e ■ A(x)G{x) dx 



± 



and so J R „ +1 V<^ • AVw = 0; thus, div A\7u = in R™ +1 . 

We now show that u is constant in all of R n+1 . Fix some X, X' G R™ +1 . By the 
De Giorgi-Nash-Moser estimate (|2.16l) and by the Poincare inequality, if r is large 
enough then 



\u(X)-u{X')\<Cr( |X X>i ) [j- |V U |2 



1/2 



By definition of N± (Vu) we have that 

\u(X) - u(X')\ < Cr 1 -^ 2 ( ^ ~ ) " (/ ^V+(Vu) 2 +iV_(V U ) 2 y /2 

and so, taking the limit as r — > oo, we have that u is constant in R™+ , as desired. 

□ 
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3.1. The fundamental solution. We now discuss the second-order fundamental 
solution. Let 2* = 2(n + l)/(n - 1), and let r 1,2 (K n+1 ) be the space of functions 
u G L 2 *(R n+1 ) that have weak derivatives Vu that lie in L 2 (M. n+1 ). From [37], we 
have the following theorems (essentially their Theorems 3.1 and 3.2). 

Theorem 3.12. Assume that A and A* are elliptic and satisfy the De Giorgi- 
Nash-Moser condition. Assume that n + 1 > 3. 

Then there is a unique fundamental solution Ty with the following properties. 

• v(X,Y) = T Y (X) is continuous in {(X,Y) G R" +1 x E n+1 : X ^ Y}. 

• v(Y) = Ty(X) is locally mtegrable in R n+1 for any fixed X G R" +1 . 

• For all smooth, compactly supported functions f defined in R I1+1 , the func- 
tion u given by 



u(X) := / T Y (X)f(Y)dY 
belongs to y 1,2 (R I1+1 ) and satisfies — divAVu = f in the sense that 

AVu ■ V<p = / ftp 



for all tp smooth and compactly supported in R™ +1 . 
Theorem 3.13. T has the property 

(3.14) / AVTy ■ Vtp = ip{Y) 

for all Y G R n+1 and all tp smooth and compactly supported in R™ +1 . 
Furthermore, T A satisfies the following estimates: 

(3- 15 ) l r r(*)l < {x Z^pi 

(3.16) ||Vr^|| L2(K „ +1VB(Y , r)) < Cr^' 2 , 

(3.17) l|vr^|| LP(B(y , r)) < cr-"+(« +1 )/f if i< p <!i±l, 

for some C depending only on the dimension n + 1, the constants X, A in (|2.5[) . 
and the constants H , a in the De Giorgi-N ash-Moser bounds. 

Finally, if f G L2(n+l)/(n+3) (W l+1 ) H if oc (R rl+1 ) for some p > (n + l)/2, then 



u(X)= / T Y (X)f(Y)dY 
is continuous, lies in Y" ' (K n+ ), and satisfies 

AWu Vip= fcp 



for all tp smooth and compactly supported in 



pn+l 



We will need some additional properties of the fundamental solution. First, by 
uniqueness of the fundamental solution, we have that if A is t-independent then 
s j(x, t) = s+r \{x, t + r) for any x, y G R™ and any r, s, t G R. In particular, 
this implies that 

(3-18) d*T£ y>s) (x,t) = (-l) k d^ s) (x,t). 
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Second, by the Caccioppoli inequality and the De Giorgi-Nash-Moser estimates, 
we have the bounds 

(3-i9) \d k n+1 r$(x)\ < C "_^ 



]/■ rl ( \'i ;)/-' r^\-', C k \X - X \ a 

\X — Y\ n+k+a ~ 1 



(3.20) |d* +1 r£(x)-^ +1 r£(x')|< 



for any k > and any X, X', Y e M™ with \X -X'\<±\X-Y\. 

Finally, it is straightforward to show that if X, Y € M n+1 with X ^=Y, then 

(3.2i) ri(Y) = rf(x). 

3.2. Layer potentials. Recall that if A and A* satisfy the De Giorgi-Nash-Moser 
condition, then the double and single layer potentials are given by the formulas 



V A f(X) = - / f(y) e ■ A 1 (y)VT^ (y, 0) dy, 

JR™ 

S A g(X)= f g{y)vf{y,Q)dy. 



Notice that by (j3~2"Tj) . dwAVV A f = and div AVd^ +1 S A g = in R' l+1 \ E™. 

We will be most concerned with the case where /, g G L 2 (R"). If A is i- 
independent, then by (|3.3j) and (|3.16|) . the integral in the definition of T> A f(X) 
converges absolutely for all X € R™ +1 \ R n and all / <E L 2 (R n ). Similarly, by 
(13.151) . if t, 7^ and if g 6 L P (R") for some 1 < p < n then the integral in the 
definition of S A g{x,t) converges absolutely. If n > 3 this implies that S A g{x 1 t) is 
well-defined for all g E L 2 (M. n ). If n = 2, so the ambient dimension n+ 1 = 3, then 
iS^g is well-defined up to an additive constant for g 6 L 2 (R 2 ). That is, if n = 2 
and if 5 G L 2 (R 2 ), or more generally if g G i p (R") for some 1 < p < n/(l — a), 
then by p. 201) . the integral 

S A g{X) - S A g(X ) - f (Vf yfi) (X) - Tf y 0) (X )) dy 

converges absolutely for all X, X Q e M™ +1 \M™. We remark that by (j3~19| . if k > 1 
then 

ff(j/)9 t fe rf yi0) (x,t)rf2/ 

converges absolutely provided g £ L p (W l ), 1 < p < oo, and if > 2 then the 
integral also converges if g g L 00 ^™). We will write 

*) = y ^ .9(2/) d*T^ 0) (x, t) 

for all such g, even if S A g does not converge absolutely. 

By (|3.19p we have a pointwise bound on d^S A g(x,t). If 1 < p < oo then for all 
integers fc > 2 we have that 

(3.22) Iflf^M)! < 

This also holds for k = 1 provided 1 < p < oo and for k = provided 1 < f> < n. 

We now establish that if (i?)^ 1 and (i?)^ 1 are well-posed, then certain properties 
of the single layer potential follow. Recall from ILemma 2.241 that under these 
conditions the square-function estimate (|2.25|l is valid. By [3j Formula (5.5)], if the 
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square-function estimate (|2.25[) is valid then it may be strengthened to the following 
estimate on the whole gradient: 

(3.23) \\\tVd t S A g\f = [ \Vd t S A g{x,t)\ 2 \t\dxdt<C\\g\\ 2 LHm . 

Suppose that A is i-independent and has bounded layer potentials, meaning that 
\\N±(VSAg)\\L 2 (R™) < C||5I!l 2 (R")- (This is Formula (|2.26p ; recall that it follows 
from (|2~25l) .) By IThcorcmTFl the operators S^'^ and (VSa)^ 1 are well-defined 
and bounded on L 2 (R"). Observe that S^' + is a Calderon-Zygmund operator by 
(|3.19|) and f|3.20|) . Thus, by standard Calderon-Zygmund theory, S^' + is bounded 
on L P (R") for any 1 < p < oo. By a standard argument (see [3TJ Proposition 4.3]), 
we may strengthen this to a nontangential bound: if g E L p (R n ) for any 1 < p < oo, 
then 

(3-24) ||JV±(^+i5Aff)||LP(R») < CWWgWupny 

The following formulas come from [3]. Again suppose that A has bounded layer 
potentials. If g e L 2 (R"), then by the proof of Lemma 4.18] we have that 

(3.25) e • A{VS A )+g - e • A(VS A )-g = -g, 

(3.26) V||5+ 5 -V||5 a5 = 0. 

In particular = regarded as operators L 2 (R n ) i-> Wj 2 (R n ). 
We now establish invertibility of and 5^ ,:t . 

Theorem 3.27. Le£ A and A* be t-independent and satisfy the De Giorgi-Nash- 
Moser condition. Suppose that A has bounded layer potentials in the sense that 
jUgl is valid. 

If (R)2 is well-posed in R™ +1 and R™ +1 , then the operator is invertible 
L 2 (R n ) 1 y W?(R n ). 

If in addition (N^)^ is well-posed in R± +1 , then S^'^ is invertible on L 2 (M. n ). 

Proof. The proof exploits extensively the jump relations for the single layer poten- 
tial. By (j3~2%1) . (l2~26|) and ITheoremTfl if g € L 2 {R n ) then 

I|3||l2(r«) = 1 1 e - A(VS A ) + g - e • A(WS A )~ g\\ 

< C\\N + (VS A g)\\ L 2 (Rn) +C\\N^S A g)\\ L 2 (WLn) . 

But if (R)% is well-posed in RJ +1 , then \\N ± (VS A g)\\ L 2^ n) < C\\ ^\\S^g\\ L ^ n) . 
By f|3.26[) we have that V\\S\~g = V^S^g and so 

We need only show that S\~ is surjective. Choose some / <E W 2 (W l ). Let u± 
be the solutions to (R)f with boundary data / in RJ +1 . By ITheorem 3.71 the 
functions g± = e • AVu±\ gR n+i exist and lie in L 2 (R n ). 

Now, let v = S A {g+ - gJ), and let v± = u| R »+i. By (f2T26|) . JV±(Vu±) <= L 2 (R"). 

Consider w± = u± + v± . By definition of u± and by the continuity relation f|3 . 26[) , 
we have that 

V||W+| aR „ + 1 =V\\W-\ gR n+i. 
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Consider the conormal derivative. We have that 

e- AWw + \ m „ +1 = e- AVu+\ m „ +1 + e- AVv + \ dR n +1 = g+ + e ■ A(WS A ) + (g+ - 9-)- 
But by (pT25l) . 

e ■ AVw+\ mn+1 =g+ + e- A(VS A )~(g+ - g-) - (g+ - g-) 
= e-AWu-\ m n+i +e-AVv_\ mn _ +1 
= e- AVw-\ ml+ i. 

Thus bv lLemma 3. Ill we have that w± is constant in RIJ_ +1 . In particular, S^(g + — 
g-) = f in W?(R n ), as desired. 

If in addition {N- 1 -)^ is well-posed in RJ +1 , then 

||fl||i=(R») < C||V||5^g|| L 2 (R „ ) < C\\N±(\7S A g)\\m«") < C||5^ ± 5||l2 (r „ ) 

and so we need only show that S^'^ is surjective on L 2 (R"). 

If (N- 1 )^ is well-posed in R!^ 1 , then for each g G L 2 (R n ) there exists some 



u with N±(Vu) G L 2 (R n ) and with d n+1 u = g on <9R± +1 . Let F = u| 9R »+i; by 
ITheorem 3.71 F exists and F G W^(R n ). Because is invertible, we have that 
F = S^f for some / e L 2 (R n ); by uniqueness of regularity solutions, S A ' f = 
g. ' " □ 

We will need some bounds on the double layer potential as well. Suppose that 
(R)i and (R)f are both well-posed in R!^ +1 , so (gUJ) is valid. By 3, Corol- 
lary 4.28], if t > then for all g £ C£°(R") we have that 

(3.28) V A (S+g)(x,±t) = T S A (e ■ A(VS A )^g)(x, ±t). 

We have the following consequences of (|3~2"g|) . Bv lTheorem 3.271 S% : L 2 (R n ) i-> 
^(E™) is invertible. Let / G Wf(R n ), and let g = (S^)" 1 /- By (|2~26t and (|3~23l) . 

(3.29) HiViCV^/)!!^^) < c\\vf\\ L2m , 

(3.30) \\\tVd t V A f\\\ 2 = ( \Vd t V A f(x,t)\ 2 \t\dxdt<C\\Vf\\ 2 L2(Rn) . 

Combining (|3 . 28f) with the jump relations (|3.25j) and (|3.26j) . we have that 

(3.31) e • AVV A f\ dK+1 - e • A\7V A f\ 9Rn+1 = 0, 
(3-32) V||2?+/-V||2? A / = -/ 
for all / G W?(R n ). 



4. The Caccioppoli inequality and related results 

The Caccioppoli inequality for second-order elliptic equations (|Lcmma 3.1[) is 
well known. In this section, we will prove a similar inequality for weak solutions to 
the fourth-order equation L* B {aL A u) — 0. We will also prove fourth-order analogs 
to some other basic results of the second-order theory. 
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Theorem 4.1. Suppose that L* B (a Lau) = inB(X,2r) in the sense oj \Definition 2.3[ 
Suppose that A and B are elliptic in the sense of (|2.5|) and that a is accretive in 
the sense of (|2.6[) . Then 



JB(X,r) r JB(X,r) r J B(X,2r) 

where C depends only on the constants X, A in (|2.5j) and (|2.6j) . 



\u\ 2 



Proof. Let ip be a real smooth cutoff function, so that ip — 1 on B(X,r), cp is 
supported in B(X,2r) 7 and |Vy| < C/r. 
Then, for any constant ci, we have that 



/" </? 2 |Vw| 2 <jRc J tp 2 Vu ■ AVu 



— Re /" <p 2 u L^u — Re / 2^^ V<^ • AVu 

and so 



^| Vu | 2 < / h 2 + L / ^ |Lau|2 _ 

J B (xar) ci 7 

Now, recall that aLAU — v for some v € Wp j oc (I?(-X", 2r)). Therefore, 

(4.2) |L A u| 2 < L~a~uL a u = j Re(v L A u), 

(4.3) H 2 <^Re(^,)<^Re(,L^). 

Because ^ is compactly supported, the weak definition of Lau implies that 
Jtp 4 vL A u = J V(<p*v) ■ AVu = J 4tp 3 vVip ■ AVu + J tp* Vv ■ AVu 



and so for any C2 > we have that 

\v\ 2 +C J ^ 2 |V(^| 2 |Vm| 



/ <p 4 vL A U < ^2 / r : J '^ 



+ 2^/^ |Vw|2 + 7?/ V ' 2|Vw|2 
<iR e y ^L X «+^ I ^ 6 |Vz;| 2 + C(1 + C2) | ^ 2 |V U | : 

This implies that 

Re J <p 4 vL A u < C(1 + C2) | ^ 2 |V U | 2 + ^ ^ 6 |V«| 2 . 
But div_B*Vt> = in the weak sense, so 

= / V(</u) • B*Vv. 
As in the proof of the second-order Caccioppoli inequality, and by (|4.3p , 
/ ^ 6 |Vv| 2 < C J ^ 4 |V^| 2 |v| 2 < CRe J Lp A \\7p\ 2 vL A u. 
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So 



Re / ip v Lau < 



J tp 2 \Vu\ 2 + ^Re J <p 4 vL A u. 



Choosing C2 large enough, we see that 

n r 



j V 4 \L A u\ 2 < jReJip 4 vL A u< ^ J tp 2 \Vu\ 



and so 



f^\Vu\ 2 <^plf \u\ 2 + ^ U 2 \Vu 

J 2r JB(X.2r) C l J 

Choosing c\ large enough lets us conclude that 

/ ^\Lau\ 2 < §/^|V«| a < £ J 



>B(X,2r 



and so 

f M 2 + l/ \Vu\ 2 <^[ \u\ 2 

J B(X,r) r JB(X,r) r JB(X.2r) 

as desired. □ 

We now prove Holder continuity of solutions under the assumption that A, A* 
and B* satisfy the De Giorgi-Nash-Moser condition. We begin with the following 
De Giorgi-Nash estimate for solutions to inhomogeneous second-order problems. 
This estimate is well-known in the case of real coefficients; see, for example, [3"2l 
Theorem 8.24]. Given the fundamental solution of [37] it is straightforward to 
generalize to complex coefficients. 

Theorem 4.4. Suppose that A, A* satisfy the De Giorgi-Nash-Moser condition. 
Let (n + l)/2 < p < 00 and let j3 = min(a, 2 — (n + \)/p). 

Then ifdivAVu = f in B(X ,2r) for some f G LP(B(X a ,2r)), then 

\u(X) - u(X')\ < C lX ~p\ 2 -^^\\f\\ LP{B{x ^ r)) 

+c£^r(f mi 2 ) 1 " 

for allX, X' G B(X ,r). 
Proof. Let 

v(X)= [ T^(X)f(Y)dY. 

JB(X ,2r) 

Let w(X) = u{X) - v(X). Bv lThcorcm 3.131 div AWw = in B(X ,2r) and so 
because A satisfies the De Giorgi-Nash-Moser condition, 

1/2 

|2 \ 



I v X'\ a 

\w(X)-w(X')\<C l - 1 



\w\ 

B(X ,2r) 



r \JB(X ,2r) J 
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By (f3~15]) . if X € R" then 



\v{X)\ < / |r£(X)| dY < Cr 2 -(" +1 )^||/|| L!)(s(Xoi2r) , 

and by applying (l3~15l) in B(X , 2r) n 2|X - X'|) and (l3~20l) in fl(X , 2r) \ 
B{X, 2\X - X'\), we see that if X, X' e R n+1 then 

\v(X) - v(X')\ < C ( |X ;5'^ M " +1)/P + I* ~ ^'| 2 - ( " +1)/p ) 
Summing these two bounds completes the proof. □ 
The following corollary follows immediately from Theorems 14.11 and 14.41 



Corollary 4.5. Suppose that L* B (a Lau) = in B(Xo, 2r) in the sense o f\Defi,nition 2.9\ 
where a, A, B are as in \Theorem and where A, A* and B* satisfy the De 
Giorgi-Nash-Moser condition. Then 

(\ !/ 2 
/ H 2 ) 
JB(X ,2r) J 

provided X, X' £ B(X ,r). 

We now prove the following pointwise estimate for solutions u in terms of their L 1 
norms. (The bound (|2 . 1 T[) is essentially the same estimate in terms of the L? norm.) 
This estimate is known for solutions to second-order equations (see, for example, 
[35| Theorem 4.1]) and may be proven for solutions to higher order equations using 
the same techniques. 

Corollary 4.7. Suppose that u is as in \Corollary 4-5 Then 



sup \u\ < C -r \u\. 

B(X a ,r) JB{X ,2r) 

Proof. Let f(p) = ||u||i«» (B(Xo(P) ) for < p < 2r. By |Corollaryl~5| if < p < p' < 
2r and if X G B(X , p), then 



B(X ,2r) 



and so 



Hp) < If(p') + ■ ( 



2 Jyr > (p> - p)«+l J BiX0}2r) ' r 

We eliminate the f(p') as follows. Let po = r and let pk+i = Pk + — r)r fe , 
where 1/2 < r n+1 < 1. Then by induction 

k — 1 

1 %. \ C f 

sup |u| = f(po) < ¥ f(pk)+^2 j- ]n+1( f \u\. 

B(X ,r) Z j=Q ( i -~~ T ) T \ ZT T Y JB(X ,2r) 

Observe that lim^oo pk = |r and so f{pk) is bounded uniformly in k. Thus, we 
may take the limit as k —¥ oo; this completes the proof. □ 

We conclude this section with the higher order analogue of (|3.3j) and some similar 
results from [3], that is, with Caccioppoli-type inequalities valid in horizontal slices. 
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Lemma 4.8. Suppose that u, <9 n +iw, and d 2 l+1 u satisfy the Caccioppoli inequality 
in that is, that whenever B(X,2r) C M™ +1 we have that 

i \Vd k n+1 u\ 2 \dt + M 2 fork = 0,1, 2. 

JB(X,r) r JB(X,2r) 

IfN+(Vu) G L 2 (R n ) then 

(4.9) BU P ||Vtt(-,t)|U a( R») <C||jV + (V«)||ia(R»). 
t>0 

// 1 > then 

(4.10) llv«(-,t)lU a (R»)<y/ H«(-,«)lli>(R»)d8 

1 Jt/2 

provided the right-hand side is finite. 
Finally, if < s < t < 2s, ffcen 

(4.11) ||V«( • ,t) - V«(« , s)\\ L2(Rn) < C— \\Vu{ ■ , r)f LHKn) drj 

provided the right-hand side is finite. 
Proof. First, we have that 

rt+KQ)/4 
Vu(x,t) — -f S7u(x,s)ds 



\Vu{x,t)Ydx < C 



dx 



/ / \Vu(x,s)\ 2 dsdx. 

Jo it 



KQ) 



But we may bound \S7u{x,t) -\7u(x,s)\ by j\Vd r u(x, r) dr\. Applying the Cac- 
cioppoli inequality to Vd r u(x,r) yields that 

(4.12) / \Vu(x,t)\ 2 dx < — — / / \Vu(x,t)\ 2 dtdx. 

Jq l \Q) J(3/2)Q it-l(Q)/3 

Applying the Caccioppoli inequality to Vu, this yields that 

r c r r t+l{Q)/2 

(4.13) / \Vu(x,t)\ 2 dx< — 3 / / \u(x,t)\ 2 dtdx. 

iQ 'WJ J2Q Jt-l{Q)/2 



Note that (|4.13[) is valid with u replaced by d n +\u. 

Dividing K" into cubes of side-length t/2, we have that (|4. 10[) follows from (|4.13|) . 
Suppose that t > 0. By definition of N+, 

f2t r 

\Vu(x,s)\ 2 dxds<Ct N + (Wu)(x) 2 dx 

lt/2 is,™ 

and applying ()4.12|) in cubes of side- length t completes the proof of (|4.9j) . 
Finally, let < s < t < 2s and let Q be a cube of side- length t — s. Then 



\Vu(x,t) - Vu{x, s)\ 2 dx 



Vd r u(x, r) dr 



dx 



<(t-s) / \Vd r u(x,r)\ 2 drdx. 
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We have that < t - s < s. Applying (|4.13|) yields that 




\d r u(x, r)\ 2 dr dx 



and so (|4.11l) is valid. 



□ 



5. The potentials £h and Jh 



We will construct solutions to the fourth-order Dirichlet problem (|2.12|) using 
layer potentials. Specifically, our solution u will be given by u = —T>aJ — Sa9 + 
^B,a,Ah for appropriate functions /, g and h. The behavior of the second-order 
potentials DaJ an d 5a<? is by now well understood fsee lSection 31 or the extensive 
literature on the subject). It remains to investigate £ = £b.u.a- Observe that 
by the definition (|2.28p of £s.a.A, when convenient we may instead investigate the 
potential T = J-B.a.A- 

In this section, we will show that J-h (and thus £h) are well-defined in R" +1 
for appropriate h and will establish a few useful preliminary bounds on J-h. Wc 
will also investigate the behavior of £ across the boundary; that is, we will prove 
analogues to the jump relations (|3.25|) and (|3.26j) . In Sections |6] and [7j we will 
establish somewhat more delicate bounds on J- (and £); specifically, our goal in 

these three sections is to show that ||iV±(V£ /i)IU 2 (R") — C||^IIl 2 (R™)- 

In lScction 81 we will show that the map h H> d n +i£h\ gRn +i is invertible L 2 (R") H> 

L 2 (W l ). We will need the assumptions that a is real- valued and A — B, or that 
||Ima||Loc and \\A — B\\l^ are small, only in ISection 81 the bounds of Sections 
[6] and [7] require only that a be accretive and that A and B* satisfy the single layer 
potential requirements of IDcfinition 2.271 We will conclude this paper by using 
these boundedness and invertibility results to prove existence of solutions to the 
fourth-order Dirichlet problem. 

We begin by establishing conditions under which J-h exists. Like the single layer 
potential Sa9, in dimensions n + 1 > 4, the integral in the definition of Fh{x,t) 
converges absolutely whenever t ^ and h G L 2 {W n )\ in dimension n + 1 = 3, Fh 
is only well-defined up to an additive constant if h £ L 2 (M. 2 ). 

More precisely, we have the following. 

Lemma 5.1. Suppose that a, A and B are t -independent, a is accretive, A and A* 
satisfy the De Giorgi-Nash-Moser condition, and Sb* satisfies the square-function 
estimate (|3.24p . 

If h G L P (M") for some 1 < p < n, and ift ^ 0, then the integral in the definition 
(|2.29|) of J^h{x, t) converges absolutely. Furthermore, 



If h € L p (W l ) for some n < p < n/(l — a), then Th is well-defined in R™ +1 and 
K" +1 up to an additive constant; that is, if we write 



(5.2) 



\Th( x ,t)\ < citl 1 -"/^!!^). 



Fh(x,t) - Fh{x',t') 




1 



d s S B *h(y,s)dsdy 
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then the right-hand integral converges absolutely. If (x',t') G B((x,t),\t\/4:), or if 
(x',t') — (x,—t), then 

(5.3) |^(a; J t)-^/i(a:' J t , )l<C'(p)l*| 1 ~ n/, 'll'»IUp(R»)- 
Proof. Recall that 

J r h(x,t)=- [ [ d s Tf s) (x,t) -j-r-d s S B *h(y,s)dsdy. 
Jr» Jo ( ' a{y) 

Choose some (x,t) G R n+1 with t ^ 0. Let B = B((x,t), \t\/2), A = A(sg, \t\/2). 
By dgUl) , \d t T± s) (x,t)\ < C/\(x,t) - (y,s)\ n . If t > and the aperture a in the 
definition of nontangential maximal function is large enough, and if 1 < p < oo, 
then 

( 5 ' 4 ) / U * ° ( \\n \dsS B 'h(y,s)\dady 
Jb KM) - (y,s)\ n 

< Ctj N + (d n+1 S B ,h)(y)dy < Ct 1 -"/ p ||iV + (a il+1 5 B ^)|| LP(R „ ) . 
Observe that 

(5 - 5) L + ^ B \( X ,t)-( y ^r^ dsM{y ' s)ldsdy 

< / N + (d n+1 S B ,h)(y) / — — dsdy. 

Jr- Jo \x - y\ m + (s + \t\) m 

If m = n and 1 < p < n, then these integrals converge and are at most 

By (|3.24p we may bound 1 1 N + {d n+ \S B * h)\\ LP r Rn -j and (15. 2[) is proven. 

To establish (JOJ), recall that by (1X1^) and ([3"^D]) , if - < |i|/4, or 

if t > and (a;', t') = (x, -t), then for all (y, s) € R" +1 \ B we have that 

Letting m = n + a, we see that if 1 < p < n/(l — a) then the right-hand side of 
(f575|) converges and is at most C(p)|t| 1 -"/ p ||^ + (a n+ i5 B ./i)|| LP(Rn) . By (pT24l) . and 
since (|5.4I) is valid for all 1 < p < oo, this completes the proof of (|5.3p . □ 



Next, we consider VFh and LaJ^H. 

Lemma 5.6. Let a, A and B be as in \Lemma 5.1[ If h G L P (M") for some 
1 <p<n/(l — a), t/ien V-7-Ti G if oc (]R!( : +1 ), awrf in particular, div AVFh is a 
well-defined element ofW^ 1 ioc (R± +1 ). 
Furthermore, 

(5.7) -adwAVFh = d 2 n+l S B ,h ml^ 1 , 

(5.8) divAVJTi^O inIT +1 
in the weak sense. 

As an immediate corollary, Eh is well-defined and also lies in W r 1 2 ;oc (RJ +1 ), and 
divAVSh = div AVFh in MJ +1 . 
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Proof. Fix some (xo,to) € R n+1 with to ^ 0, and let B r = B((xo,to),r). Let i] be 
a smooth cutoff function, supported in Bu \m and identically equal to 1 in B^u, 
with < r) < 1, |Vr?| < C/ 1 to | • For all (x,t) <E Bu i/g, we have that by definition of 
Th and by (l3~18j) . 



Th(x, t) — Th(xo, to) 

T tv s) ( x > *) TT 5 « (viV, s)d s S B *h(y, s)) ds dy 

- / r (t gi^o^o) -7-r9 s (?7(y, s)d s S B *h(y, s)) dsdy 



t (a s rf ys) (as, t) - d s rf ys) (x ,t )) 1 ^' s) d s s B , s ) d s dy 



/(x,t)- I(x ,t ) + II(x,t). 



li to < then 7 = 0. Otherwise, by (|3.22[) . the function d n+ i(r] d n+ iS B *h) is 
bounded and compactly supported, and so by IThcorcm 3.13[ I € F 1,2 (R n+1 ) C 
VK 1 2 ioc (R™+ 1 ). Furthermore, -div AVI = (l/a)9 n+ i(r? d n+1 S B »h), and so 

-a div AVI = dl +1 S B *h in B\ to \ /8 . 

We must show that V// € L 2 (-B| t[) |/ 8 ) and that div^4Vi7 = in £|t |/8- By 
(|3.20[) and lLcmma 3.11 we have that if (y, s) ^ -B|t |/4 then 

Thus by ([531) . if 1 < p < n/(l - a) then 



B |tol/ s 



I V x , t 5 s rf s) (a, t) I 1 |S S< S B * h(y, s) \ ds dy dx dt 

n+i ' \a(y)\ 

<C(p)\t \ n - n ^\\N + (d n+1 S B *h)\\ LP{ 



Thus by Fubini's theorem, 

(5.9) vufo i) = / v s , t 0,r£ s) (x, t) 1 ~^ |S) a s 5 B * %, «) ds d y . 
Jr" +i a \y) 

If <p is a test function supported in Si io i/g, then again by (|3.20[) . ()5.5[) and the 
Caccioppoli inequality, 

ip{x, t) V II{x, t) dx dt 

<C|| V || £a(R - + i ) |to| 1/2+n/a - n/, '||JV+(^+i5B./i)|| w( R») 

and so VII{x,t) G £ 2 (S| to | /8 ). 

Finally, by the weak definition (|2.8I) of div and by the formula t|5 ,9|) for V//, 
we have that div^4Vi7 = in _Bi to |/ 8 , as desired. □ 

We conclude this section by proving the continuity of V£h across the boundary. 
This property is analogous to the continuity relations (|3.26[) for the single layer 
potential, and is the reason we will eventually prefer the operator £ to T . 
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Lemma 5.10. Suppose that a, A and B are t- independent, a is accretive, A satis- 
fies the square-function estimate (|2.26j) . and B* satisfies the single layer potential 
requirements of \Definition 2.27\ Then there is a dense subset S C L 2 (R") such that 
if h G S, then 

lim \\VSh(- ,t) - V£h(- ,-t)\\ L 2 (Rn) = 0. 

Proof. We define the set S C L 2 (R n ) as follows. By assumption, Sg', + is invertible. 
Let S = (Sgl + )~ 1 (S I ), where h G S' if h[x) — u{x,t) for some t > and some u 
with divB*Vu = in R'j +1 and N+u G L 2 (R"). 

We first show that S is dense. It suffices to show that S' is dense. Choose some 
/ G L 2 (R n ). By well-posedness of (D)§ , there is some u with divi?*Vu = in 
M™ +1 and u = f on <9R™ +1 . Define f k G L 2 (R") by / fc (a;) = 1/fc); then / fc G S' . 
Bv lThcorcmT7l / fe / in L 2 (R"), and so S' is dense in L 2 (R"). 

Suppose that h G S. Then Sg« + h(x) = u(x,t) for some r > and some 
solution u. Let v(x,t) — dtu(x,t + t). By (|3.24[) . and by uniqueness of solutions 
to (-D)f , we have that dtSh = v in R™ +1 . But by the Caccioppoli inequality 
and the De Giorgi-Nash-Moser estimates, we have that N + {d n +\v){x) < ^rN + u(x) 
(possibly at the cost of increasing the apertures of the nontangential cones). 

So if h G S, then 

N + {d 2 n+1 S B *h) G L 2 (R"). 
Recall that by (|2~2"g|) and ILcmmalTTl if h G L 2 (R") and n > 3 then 

£h(x,t)=- / d s Tf Ax,t)— ^-d s S B *h(y,s)dyds- Sa( -S^l + h )(x,t). 

If n = 2 then we must instead work with £h(x, t) — £h(x, —t). 

If h G S, so that N + (d 2 +1 S B *h) G L 2 (R"), then we may integrate by parts in 
the region < s < r for any fixed r. Observe that the boundary term at s = 
precisely cancels the term S A ((l/a)S B l + h). We conclude that 

£h(x,t)= [ I Tf ) (x,t)-^—d 2 S B -h(y,s)dyds 

Jo Js. n a \v) 

r ty,r)( x ^)-^ d n+iS B 'h(y,T) dy 



f 1 

J d s Tf ys) (x,t)-j-^d s S B >h(y,s)dyds. 

Applying f|3 . 18[) and the definition of Sa , we see that 

£h(x,t) — [ Sa( —d 2 S B *h(s)) (x,t — s) ds — Sa( —d T S B 'h(T))(x,t — t) 
Jo \a J \ a J 

-J d s S A (^d s S B *h(s)^(x,t ~ s)ds. 

Here we have adopted the notation that d s S B *h(s)(y) = d s S B *h(y, s). 

Thus V£h(x,t) — V£h(x, —t) has three terms. Suppose that s £ I, s / 0. By 
\\d n+1 S B *h{s)\\ L 2 (Rn) < C\\h\\ L 2 (Rn) . Suppose that g G L 2 (R"). By gD 
and by (|2"^5)l . 

\\VS A g(s)\\ L 2 (Rn) < C\\N ± (VS A g)\\ LHRn) < C\\g\\ L 2 {Rn) . 
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Furthermore, if < t < |s|/2, then by (|4TT|) . (TO]) and ([2^26]) . 

\\VS A g{s +t)~ VS A g(s - t)||i2 (K „) < C — ||g|| L 2( R „). 

Finally by (|4.11[) . (14. 9[) . (j2.26[) and applying the Caccioppoli inequality to d n+ iS A , 
\\Vd s S A g(s + i) - Vd s S A g(s - i)IU=(R») < Cj^ja IWI^CR")- 
Thus, if t < t, then 

:>2 



\\V£h(t) - V£h{-t)\\ L 2 {m < C / \\N + (diS B ,h)\\ L 2 (wtn) ds 



C-\\h\\ m 



f 00 t 

c / — II^IIl2(R") ds. 

J T S 



Choosing r = ^/t and recalling that ft € S, we see that the right-hand side goes to 
zero as t — > + , as desired. □ 

6. A SQUARE-FUNCTION BOUND 

Recall that we intend to construct solutions to the Dirichlet problem (j2 . 12[) by 
letting u = u + + £h for some appropriately chosen function h £ L 2 (W l ). To 
prove Theorems 12 . 191 and 12 . 201 we must have that the norms || N + (W £h)\\ and 
|p V<3t£/i||[+ arc appropriately bounded. 

In this section, we will prove a preliminary square-function estimate; we will 
prove the full estimate || A r + (V£ft)||i2( R r 1 - ) + \\\t Vdt£ h\\\ + < C||ft||L 2 (R") in the next 
section. 

Theorem 6.1. Suppose that a, A and B are t -independent, a is accretive, A, A* , 
B and B* satisfy the De Giorgi-Nash-Moser condition, and S A , Sb* satisfy the 
square-function estimate (|2.25D . 

Then for all h € L 2 (M. n ), we have the bound 



\d 2 t Fh{x,t)\ 2 \t\dxdt < C\\h\\ 2 LH 



In the remainder of this section, let u = d n +iSs*h. Observe that 

d 2 Fh{x, t) = d 2 [ c> t rf " (y, s) -L u (y, s) dy ds. 

By (|2~2^| and (|X2^|) it suffices to prove that 

(6.2) \\\td?ni\\\ ± < C\\N + u\\ L 2 (Rn) + C|||tVu|||+. 

This theorem is the technical core of the paper. The proof is inspired by the T(l) 
theorems of [22] and [57]. We may think of d 2 J 7 h(x,t) as Tu(x,t) for a singular 
integral operator T with kernel (l/a(j/)) dfT^ x t \(y, s). The T(l) theorem of Semmes 
was proven by analyzing Tu — T(l)P t u and T{l)P t u for an averaging operator P t . 
We will use the same argument; our operator P t will be averages over dyadic cubes. 
Our bound on T(l)P t u will follow from well-known Carleson-measure properties of 
solutions to second-order equations. To bound Tu — T(l)Ptu in terms of Vu, we 
develop an argument ultimately allowing us to exploit the Poincare inequality. It 
bears some resemblance to the arguments of [3] (see, in particular, Lemma 3.5(i)), 
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but the particular singular integral operator at hand is different form those in [3], 
and new ideas are required. 

For ease of notation we will prove (|6.2j) only in the upper half-space (that is, 
only for \\\t <9 t 2 J-7i||| + ); the argument in the lower half-space is similar and simpler. 

Let W be the grid of dyadic Whitney cubes in R" +1 . That is, 

W = {Q = Qx [£(Q), 2£(Q)) : Q C 1" is a dyadic cube}. 

Then R™ +1 = Uq^wQ, and any two distinct cubes Q, R £ W have disjoint interiors. 
For any Q £ W we let F(Q) = {R £ W : dist(i?,Q) > 0} be the set of cubes a 
positive distance from Q and let N(Q) = Ur£W\f(Q)R be the union of cubes 
adjacent to Q. 
Observe that 



\td?Fh\\\+ = E / d t f n+1 d ^(L)(y, s) u(y, s) dyds 



If R £ W, let un — j R u. Define the four quantities 



t dx dt. 



E / ( E / d t v {I,t)(y> s )—r-;{ u (y> s )- u R) dydsj tdxdt, 



QeW"* ReF[Q) 



"=E L( E ^ 

111 = E / d t I 9 * r £) (y > s ) -tt s ) - u «) d y * 

ntliQ Jn(q) a{y) 



dy ds I i da; di, 



^= Em 



t dx dt, 



T 1 

d t T (x,t)(y^)-7-rdyds 

n(q) a (y) 



dfr^ t) (y,s)^-dyds 



i dx dt. 



We have that 



\td$Th\\\+ < C(I + II + III + IV). 



We will bound each of the terms /, II, III and IV. We begin with term IV. 

Lemma 6.3. If a and A are t -independent, A and A* satisfy the De Giorgi-Nash- 
Moser condition and Sa satisfies the square-function bound (|3.23j) . then 



IV <C\\N + u\\ 2 L2(VLn) \\l/a\\ 



Proof. Recall that 



IV = E I 2 / df f d t r£ t) (y, s) -L dy ds 
0^ JQ Jn(Q) a{y) 



1 



d i t T^ t) (y,s)—dyds 



t dx dt. 
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By the decay estimate (|3.19|) . both of the innermost integrals converge absolutely. 
By (EH, d t r£ t) (y,s) = -a a r£* t) (y,a), and so 

d 2 t I d t Tfl t) (y, s )l d yds+[ d?r£ t) (y,s)-±r dy ds 

Jn(q) a {y) Jr-+ i \jv(q) a (y) 

fl?r£ t) (y,0) -Ldy = d 2 S A {l/a){x,t). 

We claim that because 1/a £ L°°(M. n ), we have that 

d(i(x,t) = \d 2 S A (1/a) (x,t)\ 2 tdxdt 

is a Carleson measure with Carleson norm \\fi\\c at most C||1/o||£oo(r»v This 
follows from the square-function bound (|3.23[> and the decay estimate (|3. 19[) by a 
simple argument due to Fefferman and Stein (see the proof of Theorem 3 in [27]). 
Let R C W 1 be a cube; then 



i{R) 



\d 2 S A (1/a) {x,t)\ 2 tdx dt 



a 

MR) 

<2 / / \8 2 S A (l 2R (l / a)) (x,t)\ 2 tdxdt 



R 

\d 2 S A ((l-l 2li )(l/a))(x,t)\ 2 tdxdt. 

J R 



The first integral is at most C\\l/ a\\ 2 L2(2 ~ ;) < C\R\ ||l/ffl|||ao( R ») by ([5723]) . while the 

second is at most C\R\ ||l/a|| 2 oo (R „) by (|3~22|) . 

Let uw{x,t) = uq whenever (x,t) £ Q. Observe that N + uy/{x) < N + u(x) 
(possibly at a cost of increasing the aperture of nontangential cones). Then 

IV = [ \uw(x,t)\ 2 \dfS A (l/a)(x,t)\ 2 tdxdt. 

Applying duality between Carleson measures and nontangentially bounded func- 
tions, we see that 

IV < \\N + ((u W ) 2 )\\l^)\Hc < CHJV+ttllia^Hl/olliac^) 
as desired. □ 

Next, we bound the term /. 

Lemma 6.4. // a and A are t-independent, a is accretive, and A and A* satisfy 
the De Giorgi-Nash-Moser condition, then 



I < C 

Proof. Recall that 



/ |Vw(y,s)| 2 sdyds. 



1= [ ( [ 9 t T (x,t)(y> s ) -f-A u (y^ s ) - u r) dyds) tdxdt. 



QeW** r.ef(q) 
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By the Poincare inequality, 

\u{y,s)-u R \dyds<Cl{R) [ \Vu(y,s)\dyds. 

Jr 

By (|3.19p . and because £(R) < s < 2£(R) for any (y, s) £ R, we have that 



Qeff y v i?eF(Q) 



C£(i?) || l/o|| 



dist(Q,i?) 



n+2 



|Vit(y, s)| (iy <is I tdxdt 



<C\\l/a\\i 
But by Holder's inequality, 
(6.5) 



i+l (S + t + |.T - y|)"+ 2 



|Vit(y, s)| (is I tdxdt. 



+1 \7r^+ 1 (s + * + |ac — y|)' i+2 

s 2 |Vu(y, s)| 2 dyds 



V u(y , s)\ dy ds ) tdxdt 
dy ds 



+ +1 Jr^ +1 (s + t+\x-y 
C / \Wu(y,s)\ 2 sdyds. 



\n+2 



.+1 (s + t+|x-y|) 



n+2 



t dx dt 



This completes the proof. 

We may bound III similarly. 



□ 



Lemma 6.6. Suppose that a and A are as in \Lemma 6.4\ Suppose in addition 
that B* is t-independent and satisfies the De Giorgi-Nash-Moser condition, and 
divB*Vw = m R" +1 . Then 



III < C 



Proof. Recall that 



m= E / % 

Qew J Q 

Observe that 2Q c N(Q) c 5Q. Let 



| Vu(y, s)| 2 s dy ds. 



1 



d t T (x,t) (2/. s ) -jr ( u (y> s) - u Q ) ds 



t dx dt. 



w(x, t) 



T 1 

T (x,t) (V, s ) -rr (u(y, s)-u Q ) dy ds. 

n(q) a (y) 



By dS Ul), w satisfies \w{x,t)\ < C£(Q) 2 \\l/a\\ Lao(Mn) \\u - u Q \\ L ~ mQ)) . 

Bv lThcorcniTTl aL A w = (u-u Q ) in N(Q). If divB*Vu = in N(Q), then 
L* B {aLtAw) — in N{Q). Thus, we may use lTheorem 4.11 and | Corollary 4. 5| twice 
to show that 



C 

sup \dfw{x,t)\ < — -I 
0,t)eQ l W) 



l/a|U»(Rn)||u- WQ||i<»(jv(g)) 



and so 



HI < Y] [ ( 777vrllV a IU~(R")ll u - «qIU°°(jv(Q)) J icfedi. 
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Let 

n (Q) = { x '■ i x -> t) G Q f° r some t > 0} 
be the projection of Q onto M n . Notice that n(Q) is also a cube. Let U(Q) = 
6tt(Q) x {£(Q)/4,5£(Q)). Observe that if X G JV(Q), then B(X,l(Q)/4) C C/(Q). 
If div-B*Vu = in U(Q), then by the De Giorgi-Nash-Moser condition and the 
Poincare inequality, 

sup |m- u Q | <c(-f \u-uqA' <C£{Q)(-f \Vu\ 2 ) f . 

N(Q) \Ju(Q) J \JU(Q) / 

So 

///< CHl/aHioo^) V / / \Vu(y,s)\ 2 dydstdxdt. 

QeW JQJU(Q) 

Because each (y, s) G R n+1 lies in U(Q) for at most C cubes Q G this implies 
that 

m^CHl/alli^) / |V U (j/, S )| 2 S ^d S 

JR™ + 1 

as desired. □ 

Finally, we come to the term //. 

Lemma 6.7. If a is accretive, and if A and A* are t -independent and satisfy the 
De Giorgi-Nash-Moser condition, then 

2 



+1 

+ 



((S..s, II<C\ ( I %+1 {s + t+lx _ yl)n+2 \Vu(y, S )\dyds) tdxdt. 



Note that the estimate (|6.8|) together with (16.51) above imply that 
II <C I \Vu{y,s)\ 2 sdyds. 



Proof. Recall that 

dy ds t dx dt. 



11 = E / ( E / K r S)(^)^ 



To analyze the inner sum, we establish some notation. As in the proof of 
ILcmma 6.6[ if R G W is a cube, we let ir(R) be the projection of R onto R™. 
We let P(R) denote the Whitney cube directly above R, so ir(P(R)) is the dyadic 
parent of tt{R). Let 5(Q,R) = £(Q) + £(R) + dist(ir(Q),ir(R)). If G Q and 

(y, s) e R for some i? G F(Q), then 

|(y,s) - (*M)| ~ \x- y\ +max(s,t) « 5(Q,R). 
Here £7 « V if E7 < CV and V < CU . Applying (1X15]) . we see that 

To analyze the sum over R G W, we will divide VF into a "discretized cone" over 
Q and a leftover region. Let G'(Q) be the discretized cone given by 

G'(Q) = {ReW : £(R) > £(Q), £{R) > dist s tt(Q))}. 
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Let G{Q) be the larger discretized cone given by 

G(Q) = {ReW: P(R) e G'(Q).} 

Let B(Q) = G{Q) \ G'(Q) be the lower boundary of G(Q). If R is a cube, let T{R) 
be the set of cubes below R, that is, 

T(R) = {SeW: w(S) C tt(R)}. 

Observe that if R, S € -B(Q) are distinct cubes then their projections n(R) and 
ir(S) are disjoint, and so if R e -B(Q) then T(R)DG(Q) is empty. Furthermore, if S* 
is a Whitney cube, then either S € G(Q) or S € T(.R) for a unique cube i? e B(Q). 
Thus, we may write 

"<^E /( E i^- ^J^U ) 2 ^ 



We may write |ug — uq\ < \us — ur\ + \ur — uq\ in the second sum. Observe that 
if R € S(Q) = G(Q) \ G'(Q), then i? is not in G'{Q) but the cube P(R) above it 
is, and so either i(R) w or w dist(7r(.R), tt(Q)). Thus if 5 e T(R) then 
5(5, Q) « w <5(i?, Q). So we may write 



2 

tdxdt 



+ c E/i E £ i«— isg^M 



2 

tdxdt 



Qew J Q RaB(Q) seT(R) 
= C(V + VI). 

We begin by analyzing the term VI. Choose some R e B(Q). Let Lo = {R}, 
and let Lj = {S e T(i?) : £(5) = 2"^( J R)}. Then 

OO 

E i«s-«fliisi = ifliE 2 " J ' (n+1) Ei^-^i- 

SeT(P) j=l SGL 3 

Let [/,• = X)seL J l u S' — u r\: observe that Uq = 0. Then for each j > 1, 

u j = E \ us - ur \ - E i Ms ~ + E i up ( s ) ~ 

seLj seLj seLj 

The second term is equal to 2 n Uj-\. To contend with the first term, we apply the 
Poincare inequality in the set S U P(S). Then 

\us - up(S)\ < \us - U SuP{S )\ + |u P (s) - USUP(S)I 

If, 1 



= TCI / l« - U SUP(S) I + ^T7n / \u - USUP(S) 



5| Js' 1 2«|5| 



P(S) 



i r £(5) /■ 

- Tei / I u -«sup(S)I < / |Vu(t/,s)|cfyds. 

PI JSUP(S) PI JSUP(S) 
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Thus, 

Uj < 2 n Uj-! + 

seLj " v "' uo seL J 
Since Uo = 0, by induction 

v — . « ( / « « : / 

|V«| 

is 



E^/|v»i +2 -E^/iv»!. 

1 



which may be simplified to 



t=n ficr.,_ \ / J S 



k=0 SeL k 

Thus 

OO 

E \u s -u R \\S\ = \R\J2^ j{n+1) U j 

SeT(R) 3 = 1 



< |i?|^2-^+ 1 )^2"0- fe ) ^ — - / \Vu(y, S )\dyd S 
j=i k=o seL k { > Js 

oo C f 

= \R\J2i- k - nk E j(sy JJ^(y,s)\dyds 



k=0 SEL k 

< 



SeT(R)u{R} 

So 



C E t(S) f\Vu(y,a)\dyda. 



VI = E £( E m ry^ S **** 

C E / ( E I77j mn+2 E ^) / |V«(i/,*)|dj/<to) t<*rdi 

' L +1 (/«- (t +s+ il- y i)«- |vu(y ' s) 1 dy ds ) * dx dt 



< _ 

QeW J Q^ReB(Q)° W,ri> ' SeT(R)u{R} 

2 

III s \ 

< c 

7 V — + 

because 5(Q, R) w <5(Q, 5) w i + s + |a; - y| for any S e T(i2), any i? € and 
any (a;, t) £ Q and (y, s) G 5. 
Now, recall that 



v = E / ' ( E l u «-»ql Arn l gL + 2 > ) icfa ^ 



Qew-'y y ReG(Q) ; 

Observe that if R e then dist(7r(i?), tt(Q)) < C^(i?) and so S(R,Q) w ^(i?); 

thus 

^ C E /( E I u *- u «It?W) <dTdi - 

Let Gj(Q) = {i? e G(Q) : ^(i?) = 2H{Q)}. Let x Q be the midpoint of the cube 
7r(<2). There is some constant c depending only on dimension such that if R £ Gj (Q) 
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then R is contained in the cylinder C 3 {Q) = A(x Q ,c2H(Q)) x {2H{Q), V +1 l.{Q)). 
Therefore, 



\u R -u Q \ = + \ } I \u(y,s)-u Q \dyds 

R£Gj(Q) IVI ReGj(Q) 

<Cf \u{y,s) - u Cj ( Q )\dy ds + C\u C jq) - u Q \ 

If j > 1, then 

3 

\ U C 3 (Q) ~Uq\ < \u C „(Q)-U Q \+J2\ U C k (Q) -"C^iCQ)!- 

k=l 

Applying the Poincare inequality in Ck(Q) U Ck-i(Q), we see that 

C [' 

\uc k {Q) ~ u Ck _ l{ Q)\ < 9knf(n)n / \Vu{y,s)\dyds. 
Observe that Q C C (Q), and so by the Poincare inequality, 

\u Co (Q) ~ u Q \ < C£(Q) f \Vu(y,s)\ dyds. 

JCo(Q) 

Finally, f Cj{Q) \u(y, s) - u Ci (Q)\dyds < C2H{Q) f c . {Q) \Vu(y,s)\dyds. So 

E \ur-u q \ <^___ / \Vu{y,s)\dyds. 
rgG 3 (Q) k=o Z lW) Jc ^Q) 

Thus, 

1 °° 1 



oc 1 j c 

= CV w 2 y Q) N / I Vu(y, a) | da. 



But if (y,s) e Cfe(Q) and e Q, then |a; - j/| +i < C2 k £(Q) and s w 2 fc £(Q). 

Thus, 



^ C E / ( E I u «- u «I7tW) * da:cft 
" Q^^QV^ 2 fe (" +2 ^(g)" +2 ic fc (Q) 1 y 

< C / [/ t, 1—^ — rr|Vu(y,s)| dwds ) iefecft. 



This completes the proof. □ 
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7. A NONTANGENTIAL BOUND 

We have established that, if h G L 2 (R"), then under appropriate assumptions 
on a, A and B, 

(7.1) f [ \d?Th( Xl ±t)\ 2 tdtdx <C\\h\\ 2 L2(Rn) . 

jr™ Jo 

In this section we will prove that, if h € L 2 (R"), then under appropriate as- 
sumptions on a, A and B, 

|||tVft^)||| + + \\N + (\7£h)\\ LHRn) < C\\h\\ 2 L2(Rn) . 

Recall that Eh = Th - S A {\S^h); thus by flOg) and ([5~2l?]) , the difference 
— J 7 /! satishes square-function estimates and nontangential estimates, and so we 

may work with £ or J 7 , whichever is more convenient. 
This proof will require several steps. 

Lemma 7.2. Suppose that v € W 2 ; oc (R± +1 ) satisfies the conditions of \Lemma 4-8\ 
that is, that v, d n +iv and d 2 +1 v satisfy the Caccioppoli inequality in RJ +1 . 
Suppose furthermore that there is a constant C\ such that 



(7-3) \\Vv(- ,t)\\ L 2 (Rn) <d\t\ 



-l 



Then there is some constant C , depending only on the constants in the Caccioppoli 
inequality, such that 

POO P POO P 

/ / \Wv(x,±t)\ 2 tdxdt<C / \d t v(x,±t)\ 2 tdxdt + C'l 
Jo Ju n Jo il" 

We will immediately apply this lemma to dt JFh in the lower half-space. Later in 



Corollary 7.12 we will apply this lemma to dt£h in the upper half-space as well. 



The bound (|7.3I) is simpler to establish in the lower half-space because J-h(x, t) 
is defined in terms of an integral over and certain direct bounds can be 

computed if (x, t) does not lie in the integrand. 

Proof. This parallels the proof of (31 Formula (5.5)], where a similar inequality was 
proven for the single layer potential. 

Define u m (t) — J R „ | V9 t m «(j;, ±t)\ 2 dx. Observe that by assumption uo(t) < 
C 2 t~ 2 , and by the Caccioppoli inequality, 

t 3 u 1 (t)dt=f \Vd t v(x,±t)\ 2 t 3 dxdt < C f \d t v(x,±t)\ 2 tdxdt. 

JR^ +1 JR^ +1 

We wish to bound f °° t uo(t) dt. 

Suppose t > 0. By (14.11[) . uq is locally Lipschitz continuous. Thus if < e < 
S < oo, then 

S pS pS pS 

tuo{t)dt~ I tuo{S)dt— I t u' Q (s)dsdt 



C? 1 



s 



n 
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But |u (s)| < 2y / u (s) ui(s) and 



(72 _ 

t Uo (t) dt < + / s 2 \J uq(s) ui(s) ds 



C i 1 f S / \ j 1 f S 3 
2 + 2j su °(s)ds+- J 



Rearranging terms, we have that 

•S pS 

tu Q (t)dt<Cf + s 3 u 1 (s)ds. 
Taking the limit as e — > + and S — > 00, we have that 

/>OG /* />00 />00 

/ / \\7v(x,±t)\ 2 tdxdt = tu {t)dt<C 2 + t 3 Ul (t) dt 
Jo Jm™ Jo Jo 

<Cf+C f f \d t v{x,±t)\ 2 tdxdt 

Jo il" 

as desired. □ 

We now proceed to nontangential estimates in the lower half-space. 

Lemma 7.4. Suppose that a, A and B are t -independent, a is accretive, A satisfies 
the single layer potential requirements of \Definition 2.2% B and B* satisfy the De 
Giorgi-Nash-Moser condition andSs* satisfies the square-function estimate (|2.25l) . 
Then for every h € L 2 (M. n ), we have that 

(7.5) \\N-(V£h)\\ L 2 m < C||J»|| La(B ») 

and the boundary value X7£h\ gs n+i exists in the sense of (|2 . 10[) . 

Proof. By the definitions (f2~29| and (pT23|) of T and Sa, 

(7.6) Vd t J 7 h(x,-t) = - J VdtS A (^~d s S B *h(s))(x,-t- s)ds 

where as in the proof of lLcmma 5.101 we let SB*h(s)(y) = SB*h(y,s). Because A 



and B* are i-independent and have bounded layer potentials, if t > 0, then by 
(|3.24[) and the Caccioppoli inequality, || V9 t J r /i( • , — £)||l 2 (R") < Ct _1 ||/i||| 2 j R „- ) . 
Thus bv lLcmma 7.21 

/ / \Vd t Fh{x,~t)\ 2 tdxdt<C / \d%Fh(x,-t)\ 2 tdxdt + C\\h\\ 2 L2(Rn) . 

Jo JR n Jo JR n 

Bv lThcorcm 6. II this is at most C||/i||^ 2 ( Rre ). 

Since div AX7(dt Fh) = in K™ , we may apply [Theorem 3"Jjl in the lower half- 
space. Thus 

(7.7) \\N-{d t Th)\\l HVLn) < ( \Vd t Th{x,-t)\ 2 tdxdt<C\\h\\ 2 L2 ^ n) 

and the boundary value 

J 7 - 1 ' h = d n +i .F/i| aK „ + i 
exists in the sense of L 2 functions. 
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Because A satisfies the single layer potential requirements, S A ' is invertible 
L 2 (R n ) i-> L 2 (R n ), and so there is some g £ i 2 (R") with S^~g = Fur- 
thermo re, ||g ||/^( R ») < C\\h\\ L 2 {Rn} . 

By (|3.24l) . if u = d n+ iS A g, then ||iV_u||i3rR»i < oo. By (|7.7|l the same is 
true of u = dn+iJ-h, and so by uniqueness of solutions to (Dfe, we have that 
d n +iS A g = d n+1 Th in R^ +1 . By (2~25]l . and by and ITheorem 4.T1 we have 
that m = — Th satishes the conditions of ILcmma 3.91 and so S A g = Fh in 
up to an additive constant. 

In particular, by (|2.26p we have that 

\\N_(VTh)\\ L 2 {Rn) <C\\h\\ L 2 {Rn) . 

By definition of £ h and by (]2 . 26[) , (|7.5[) is valid. By ITheorem 3.71 and because 
div AW(£h) = in IT +1 , we have that \7£h\ gRn +i exists and lies in L 2 (R"). □ 



Lemma 7.8. Suppose that a, A and B are as in \Lemma 5.1\ Let h £ L 2 (M. n ). 
Then 

\\N + (d t £h)\\ L 2 (Rn) < C\\N_(d t £h)\\ L 2 (Rn) + C\\h\\ L 2. 

Proof. This is essentially the analogue to Cotlar's inequality for singular integral 
operators (see, for example, [331 (8.2.2)]) reformulated to apply to our potential £. 
We prove it using similar arguments. 

Let x* £ M. n and let (x, t) £ r y + (x*) for some t > 0, where 7+ is the nontangential 
cone given by (12.21) . There is some constant jo, depending on the aperture a of the 
nontangential cones, such that x £ A(x*,2 J0 t). Let A = A = A (a;*, 2 J ' 0+1 f), and 
for each j > 1, let Aj be the annulus A(x*, 2^°+^ +1 t)\ A(x*, 2 jo+ H). Let hj = hl Aj , 
so h = J^jLo hj ■ Then 

00 

\d n+1 £h(x,t)\ < \d n+1 £h(x,-t)\ + ^2\d n+1 £h J (x,t) - d n+ i£hj(x, -t)\. 

3=0 



Observe that \d t £h(x,-t)\ < N-{d n+x £h){x*). By (JO} and |Corollary~4~5l if 1 < 
p < n then 

\d t Fh 3 ( Xl t)\ + \d t Th 3 (x,-t)\ < ^IIM^Rn) < CV n /PM(\h\i>)(x*) 1/p 
where M denotes the Hardy-Littlewood maximal function. Recall that 

£hj(x,r) =J r hj(x,r) - Sa (^Sgi + hj^j (x,t) 

= - / ds~Tf ){x,t) -j—dsSs^hjiy, s)dsdy-S A ( -Sgl + h, } ■) (x,r). 

JRI+ 1 {y ' ' a[y) V a / 

Applying (13.22J) to S A and p. 241) to Sgl. + , we see that the same is true of dt£hj. 
We will apply this bound only for j = and j = 1; we seek a better bound for 
3 > 2. 

Suppose j > 2. Let A, = A,-i U A,- U A7+1. We divide K™ +1 into three pieces: 
I = Aj X (0, 2^), 77 = Aj x (2^'i, 00) and III = (E n \ A,-) x (0, 00). We integrate 
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by parts in s in II and 77/. Recalling the definition of Sa, we have that 
d T £hj(x,T) = J^d 2 s Tf y s) (x,t) -^-jd s S B »h : j(y,s)dsdy 

<9 s rf } (x, t) — !— d 2 S B * hj (y, s) ds dy 
ilium a iy) 

+ J A ^1^ 2^)0, t) d n+ iS B 'hj(y, 2H) dy 



1 



~d T S A \ l A .~S B l hjj (x,r). 

Thus d n j r i£hj(x,t) — d n +i£hj{x, —t) may be written as a sum of four terms. Ap- 
plying the bound (|3.24[) on Sb* , the bound (|3.22[) and Holder continuity to Sa and 
S B * , and the bounds (|3.19l) and p.20[) on T^ y 8 Jx, r), we may show that each term 

is at most C2-i a M{\h\P){x*) 1 /P. 
Thus, if 1 < p < n then 

oo 

\d t Sh(x,t)\ <N_{d n+1 Eh){x*) + Y,C2- ja M{\h\P){x*) 1 / p . 

j=o 

Choosing 1 < p < 2 and recalling that M is bounded LP i-> LP for any 1 < p < oo, 
we have that 

\\N+(d n +i£ h)\\ L 2( Rn) < \\N-(d n+1 £h)\\ L 2 (Mn) + C||/i|| L 2 (Rre) 
as desired. □ 

Lemma 7.9. Suppose that a, A and B are t-independent, that a is accretive, and 
that A and B* satisfy the single layer potential requirements of \Definition 2.2% 
Then for all h £ L 2 (R n ), 

(7.10) \\N + (V£h) 2 \\ mRn) < C\\h\\ LHRn) 
and Vf /i| aR n+i exists and satisfies 

(7.11) V^| aB » +1 =V^| att » +1 . 

Proof. First, assume that h £ S where S is as in ILemma 5.101 By Lemmas 15.101 
and l7.4[ the formula (I7.11|) is valid, and so we may define £ ± h = £h\ gm n+i. 

Let x* £ R n and let (x,t) £ -f(x*). Bv lTheorem 4.11 and |Corollary 4.7[ we have 
that 



/ \V£h\ 2 <^-f \£h~E\ 2 <%(f \£h 

J B((x,t),t/2) r J B((x,t),5t/8) r \JB((x,t),3t/A) 



for any constant E. 
If (y,s) £ then 



\£h(y, s)-E\< \£h(y, s) - £+h(y)\ + \£+h{y) - E\ 
< sN + (d n+1 £h)(y) + \£+h(y) - E\. 
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Let A = A(x*,Ct) where C is large enough that A(x,t) C A(x*,Ct). Then 
/ \V£h\ 2 < ~ ( / tN+(d n+1 £h)(y) + \£+h(y) - E\ dy) 

J B({x,t),t/2) 1 \Ja J 

- c (£ 7V+(a " +i£/i)(y)dy ) +^{fj £+h ^- E \ dx 

Choose E = j A £h(x, 0) dx and apply the Poincare inequality. We conclude that 
N+(\7£h)(x*) 2 = sup i \\J£h\ 2 

(x,t)ey(x») J B({x,t),t/2) 

< CM(N + (d t £h))(x*) 2 + CM{V\\£ + h){x* f. 

Bv lLemma 7.4l and lLemma 7.81 and by the L 2 -boundedness of the Hardy-Littlewood 
maximal operator, we have that (|7.10p is valid for all h £ S. 
We now must pass to arbitrary h £ L 2 (M. n ). Define 

N s , R F(x) = sup( ( / |F| 2 ) :(y, S )€ 1+ (x) ) 5<s<R,\y\<R). 

By ([O]) . (pQ2)) and ITheorcm 4771 for any fixed S > 0, R < oo the map h i-> 
Ns,u(V£h) is continuous on L 2 (R"). Thus because 5 is dense in L 2 (IR n ), we have 
that for any h e L 2 (R n ), 

\\Ns. R (V£h)\\ L 2 {Rn) <C\\h\\ L 2 (Rn) 

uniformly in 5, R. Letting S and R —> oo establishes (|7.10p . Recall that (|7.11j) 
is valid for all h G S. By (|7.10p and ILemma 4.8| we may extend (|7.1ip to all of 
i 2 (R"). □ 

Corollary 7.12. Suppose that a, A and B are as in \Lemma 7.9[ If h £ L 2 (M. n ), 
then 



f f \Vd t £h(x,t)\ 2 tdxdt<C\\h\\ 2 L2{ 

JO JM" 



»)• 



Proof. By ITheorcm 6.11 (|2.25[) and the L 2 -boundedness of Sg', + , we have that 



B* 

\df£h{ Xl t)\ 2 tdxdt < C\\h\\ 2 L2{Rn 



By Lemmas 14.81 and 17.91 we have that 

\\Vd t £h( - ,t)\\ L 2 iMn) < j\\N + (V£h)\\ L 2 {win) < j\\h\\ L 2 {Rn) 
and so the conclusion follows from ILemma 7.2| □ 

8. The proof of the main theorem 

In this section, we will first prove an invertibility result for the potential £ = 
£b.o,.a- We will use this invertibility result to prove existence of solutions to the 
fourth-order Dirichlet problem. Finally, we will conclude this paper by proving 
uniqueness of solutions and establishing a Fatou-type theorem. 

We remark that the invertibility argument and the construction of solutions 
in this section closely parallels the construction of solutions for the biharmonic 
Dirichlet problem of [3D] and [S3]. 

The invertibility result we will prove is the following. 
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Lemma 8.1. Suppose that a, A and B are t-independent, that a is accretive and 
that A and A* satisfy the single layer potential requirements of \Definition 2.2% 
Then there is some e > and some C > 0, depending only on the parameters listed 
m l Theorem 2.19\ such that if 

|| Ima|| £,«.(][[») + \\A - fl||i«.(H») < e 

then the mapping h > d n +i£B, a ,Ah\g^+i is invertible on L (W 1 ), and its inverse 
has norm at most C. 

Proof. Define £g a A h — <9 n +if_B,a,A^laR"+ 1 ■ We begin by showing that £ A a A is 
invertible for a real; it is for this reason that we require that both A and A* satisfy 
the single layer potential requirements. 

Choose some h G S, where S C L 2 (R") is as in ILcmma 5.101 By definition of 
the single layer potential, 

/ g(x)d t J 7 A,a,Ah{x,-t)dx= f [ d 2 S A rg(y, s + t) —}— d s S A *h(y, s) ds dy. 
Jr™ Jr" Jo a (y) 

Integrating by parts and applying the dominated convergence theorem yields that 

/ 9(x)£A. a .Af( x )dx = - lim / / d s S AT g(y, s - t) — ^- d 2 S A * f(y, s) ds dy 
Jr" ' ' t->0- J M „ J a a{y) 

= Jr- S ^ +9 ^ 'aJy) S ^' + dy ~ f un ^ £ ^,a,Ad( x ) dx - 

By density of S and by the L 2 -boundedness of £ Aa A and £^ aA , we have that for 
any h E L 2 (M"), 

/ h£X a>A h+ f h£^ a ^h= [ -±-S^ + h(y)S^h(y)dy. 

Jr™ Jr» ' ' Jr" a\y) 

Observe that £A,a,Ah(x, t) = £ A - A h(x,t). Thus, if a is real, then 



hE ^aA hda= / ~ h£ A^A hd(J 

and so 

Re f^h£i^ A hda=\j^ -L.\S^h(y)\ 2 da(Y) > ^\\S^ + h\\l, {Rn) . 

Thus because A* satisfies the single layer potential requirements, 

1 



h£ia,A h do- 



>^\\hf L2 



and so £ A a A must be one-to-one. But the adjoint (£ A a A )* must also be one-to-one, 
and so £ AaA must be invertible. Furthermore, (^ aj4 ) _1 is bounded L 2 (E n ) H> 
L 2 (R n ), as desired. 

We prove invertibility for ||Ima||ioo or ||^4 — -B||l°° small using standard analyt- 
icity arguments. Specifically, the mapping a *-> £ Aa A is analytic in the sense that 
if z H> a z is analytic D H> L°°(M. n ) for some DcC, then z t- > £ Aa A h is analytic 
D i — ^ L 2 (M. n ) for any h € L 2 (W n ). Furthermore, under our assumptions, if ao is 
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real and accretive then \\£ A a A^W^tM 71 ) is bounded for all i-independent a in a L°° 
neighborhood of ao; thus if ||a — ao||L°c is small enough, then 

II^A.o,^ - ^i,a ^IU 2 (K")^L2(R") < C\\a - a ||Loo( R „). 

Thus, if ||Ima||x,oo/ R n) is small enough then £ A a A is invertible on L 2 {R n ). Similarly, 
invertibility of £ B a A for ||^4 — _B||^oc( Rre ) small enough follows from analyticity of 
the map B H> Tf . □ 

We now use this invertibility result to prove existence of solutions to the fourth- 
order Dirichlet problem. 

Theorem 8.2. Suppose that a, A and B satisfy the conditions of \Lemma 8.1[ Let 

f G W?QSL n ) and let g £ L 2 (R n ). 

Then there exists a constant c and an h G L 2 (M. n ), with ||/i||z, 2 (R»») < C|| V/IIl 2 (r™) + 

C||fflU 2 (R"); such that 

u{X) = -V A f{X) - S A g{X) + £ B , a ,Ah(X) + c 

satisfies 

L* B (aL A u) = mM™ +1 , 

u = f ondW l + +1 , 

e-AVu^g on <9R™ +1 

where L* B (aL A u) — in the sense of \Definition 2.9] and where u — f , e ■ A'S/u = g 
in the sense that 

Hm ||u(-,*) - f\\ W 2 (Rn) + ||e • AVu(- ,t) - g\\ L 2 {wl) = 0. 

Furthermore, 

(8.3) ||jV+(Vti)|U 3(B ») + \\\tL A u\\\ + + lll*Vft«||| + < C||V/|| L 2 (R „) + C\\g\\ L 2 {Rn) . 

Proof. Let v(X) = -V A f(X) - S A g(X), so that div AVv = in R" +1 \ E™. By 
and (ggQ , and by (1330)) and (pHB")) . we have that 

(8.4) \\N ± (Vv)\\ L 2 (Rn) + \\\tVd t v\\\ < C\\Vf\\ L 2 (Rn) + C\\g\\ L 2 {Rn) . 

Bv lThcorcm~3~7l d n+1 v\ dRn +i exists and lies in L 2 (R"). Let h G L 2 (R n ) be such 
that 

5n+lfs,a,A/i| 9K n+i = — d n +iv\ gRn + l . 

Bv lLemma 8.1[ h exists and satisfies ||/i||l 2 (r™) < C | j V/||l 2 (r™) +C||g|| i 2( R ^- ) . Thus 
by ILemma 7.91 and |Corollary 7.12| 

(8.5) \\N±(V£ B ,a,Ah)\\ L , iUn) + \\\t Vd t £ B ,a,Ah\\\ < C||V/|| i2(R „) + C\\g\\ L 2 (Rn) . 
By L A £ B ^ A h = (l/a)d 2 +1 S B *h and so by (233, 

< C||fc|U 3(B ») < C||V/|U 2(R n) + C||«7|| i3{B n). 

Let w = v + £ B , a ,Ah = —T>Af — S A g + £b,h, Ah. Then w satisfies ( )8.3|) . By (|5.7[) . 
L* B (a L A w) = in R™ +1 in the sense of lDcfinition 2.9i We need only show that for 
some constant c, u = w + c has the correct boundary values. 

By our choice of h, d n+1 w = on <9R" +1 . By (glD and ([53]) . <9 n+ iw is a 
solution in R" +1 , and so 9 n+ iw = in R" +1 . Again by ((8~4|) and ([875|> . w satisfies 
the conditions of lLcmma 3.9[ and so w is constant in R™ +1 . 
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By (1X251) . (pOf)]) . (EOTl) . and (EO^ . 

lim Viiuf- , *) — V|iu(-,-t) = V/, lim e • AVv{ ■ , i) - e ■ AVv(- , -t) = g 

in L 2 (R n ), and by (|7. 1 1[) the same is true of w. Since w is constant in this 
yields that e • iVw| 3K n+i = g and that V||w| aR n+i = V/. As in IRemark 2.181 we 

have that it? — > / — c, for some constant c, pointwise nontangentially and in L 2 (M n ). 
Letting u = w + c completes the proof. □ 

We conclude this paper by proving a Fatou-type theorem and uniqueness of 
solutions. 

Theorem 8.6. Suppose that a, A and B are t -independent, that a is accretive, and 
that A and B* satisfy the single layer potential requirements o f\Definition 2.2% We 
do not require that ||Ima||ioc or \\A — -B||l°° be small. 
Suppose that u satisfies 

( L* B (aL A u) = mRl +1 , 

\ ||jV + (Vu)|| ia( R„) + |||*Imu|||+ < oo. 

Then G = Vu| aR n+i exists in the sense of (12. 10|) and satisfies 

(8.7) l|G|| W) < cuiv+^iu^Kn) + C||[ti A u||| + . 

Furthermore, there exist functions ua & W 2 ; oc (IR™ +1 ) and h € L 2 (M. n ) that 
satisfy the estimate 

\\N+(Vu A )\\ LHRn) + \\h\\ L - HRn) < C\\N+(Vu)\\ LHRn) + C|||tL A u|||+ 
and such that 

u = ua + ^B.a.Ah and div AVua — in 

Proof. If L* B (a Lau) = in the sense of IDefinition 2.9| then there is some w € 
j oc (K. ) such that w = a Lau in the weak sense. Furthermore, div£?*Vw = 
inK™ +1 . 

Now, observe that by the De Giorgi-Nash-Moser condition. ITheorem 4. H and the 
Poincarc inequality, and the definition of JV+, 

\w( x ,t)\<c -f H 2 <? \i iv^l 2 



' B((x,t),t/4) J t \J B((x,t),t/2) 

< t , (V) < ? ( / (V«) 

t |x-i/|<t/C t V J A(x,t/C) 
< Ct -l-n/ 2|| ^ +(Vu)||L2(Rji) 

We define 

i) = / w(x,s)ds= I w(x,t + s)ds 
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Observe that the integral converges absolutely for all t > 0. Furthermore, v <E 
w i.ioc(M+ +1 ) and L B'V = 0. Using ILemmaTSI we may bound Vv as follows: 

||Vu(-,t)||i,a(R») < C / ||Vw(- ,s)||i2 (R „) ds < / — \\w(- ,s)|| L 2 (R „) ds 



- / ^\\N + {Vu)\\ L 2 {Rn) ds< j\\N + (Vu)\\r* {R ny 



Thus, by ILemma 7.21 



\Vv(x,t)\ 2 tdxdt<C / \d t v(x,t)\ 2 tdxdt + C\\N + CVu)\\j 2m n). 

JRI+ 1 { ' 

But since d n +\v — ~w — —clLau, the right-hand side is at most 
C\\N + (Wu)\\ 2 L2(Rn) +C\\\tL A u\\\l. 

Thus, by ITheorem 3.61 

||JV+«|U»(*») < C\\N + (Vu)\\ L * m + C|||ti A ti|||+. 

By invertibility of S^', + and by uniqueness of solutions to (D)^ , we have that 
v = —d n +iSB*h for some h € L 2 (W l ). So clLau = w = d 2 +1 SB'h. 

Let ua = u — Ss.a.Ah. By ILemma 5.61 divAVuA — in Furthermore, by 

ILemma 7.91 

\\N + {Vu A )\\m^) < \\N + {Vu)\\ L 2 {Rn) + C\\h\\ L 2 (Rn) 

as desired. 

The existence of G = Vu\ aM n+i and the bound (|8.7p follows immediately from 
ITheorem 3.71 and ILemma 7.91 □ 

Finally, we prove uniqueness of solutions. 

Corollary 8.8. Let a, A, and B be as in I Theorem 8. 61 Assume in addition that 
^B,a,A * s one-to-one, where Sg a A h = d n+ i£B,a,Ah\ gm ^+i . 
If u satisfies 

L* B (aL A u)=0 inWl + \ 
Vu = on dR+ +1 , 
, \\N+(Vu)\\ L 2 {Rn) + ||[tiA«|||+ < oo 
then u is constant in M7 + . 

Proof. By ITheorem 8.61 we have that u — ua + £B,a,Ah for some h £ L 2 (M. n ) and 
some ua with divAVuA = in 

Let u + = ua, and let it_ = —£s,a,A n in Observe that divAVu± = in 

Bv lLcmmaTJl and because Vu = on we have that 

Vu_] „ +1 = -V£ B ,a,Ah\ mn+ i = -V£B,a,Ah\ m „+i = Vu + \ m „+i. 



Thus bv ILemma 3.1l| u+ = ua and U- = —£B,a,Ah are constant in In 
particular, if £^ a A is one-to-one, then h = 0, and so u = ua + £B,a,Ah is constant 
in as desired. □ 
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